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ABSTRACT 


Two  methods  for  evaluating  the  accuracy  of  hydrographic 
positioning  data  are  presented.  One  method  consists  of 
classifying  each  position  in  a  survey  based  on  the  radius  cf 
the  90  percent  confidence  circle.  The  second  method 
involves  classification  of  positions  based  on  the  parameters 
of  the  90  percent  confidence  ellipse.  Both  methods  are 
based  on  geometric  and  statistical  relationships  between 
intersecting  lines  of  position. 

Range-range,  azimuth-azimuth,  and  range-azimuth  posi¬ 
tioning  data  are  classified  using  both  criteria.  For 
noncritical  positions,  the  confidence  circle  method  is  found 
to  be  preferable  due  to  its  ease  of  interpretation.  For 
positions  of  significant  features,  such  as  underwater 
hazards,  the  confidence  ellipse  provides  a  more  useful 
representation  of  the  shape  and  orientation  of  the  true 
error  distribution. 

The  concept  of  presurvey  positioning  design  is  alsc 
presented.  With  the  aid  of  computer  graphic  displays,  the 
hydrographer  can  predict  the  accuracy  of  offshore  posi¬ 
tioning  data  prior  to  data  acquisition.  By  analyzing  accu¬ 
racy  lobes  generated  about  shore  stations,  a  survey  can  be 
designed  to  meet  given  specif ications. 
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I.  INTRODUCTION 


8.  BACKGROUND 

A  hydrographic  record  can  he  viewed  as  the  resultant  of 
two  independent  measurements  made  at  a  discrete  point  ever  a 
tody  of  water.  These  measurements  involve  the  determination 
of  a  vessel’s  position  at  a  given  time  as  well  as  the  depth 
of  water  at  that  position.  Of  interest  to  the  hydrographer 
and  to  the  user  of  hydrographic  data  is  the  accuracy  of  the 
position  deter minatiens.  Fundamental  to  the  determination 
of  positional  accuracy  is  the  identification  of  the  sources 
of  errors  in  position  measurements  and  the  ultimate  treat¬ 
ment  of  these  errors. 

A  hydrographic  position  can  be  determined  by  a  number  of 
methods  all  involving  geometric  relationships  between  known 
points  and  the  vessel’s  unknown  location.  The  known  points 
may  be  fixed  stations  on  shore,  whose  coordinates  have  teen 
determined  by  geodetic  survey  methods,  or  they  may  be 
rapidly  moving  satellites  whose  coordinates  in  time  and 
space  can  he  defined  very  precisely.  A  hydrographic  posi¬ 
tion  is  established  by  the  intersection  of  two  or  mere  lines 
of  position  (LOP's)  which  are  generated  by  the  geometric 
relationships  between  the  fixed  points  and  the  vessel's 
unknown  location.  The  resultant  accuracy  of  the  vessel's 
position  is  therefore,  in  part,  a  function  of  the  errors 
associated  with  the  intersecting  LOP's. 

Several  measures  cf  accuracy  can  be  used  to  evaluate  the 
quality  of  a  hydrographic  position.  Predictability,  or 
absolute  accuracy,  is  the  measure  of  accuracy  with  which  the 
positioning  system  can  define  the  location  of  the  same  point 
in  terms  of  geographic  coordinates  Repeatability,  cr 
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relative  accuracy,  is  a  measure  with  which  a  positioning 
system  permits  a  usei  to  return  to  a  specific  point  on  the 
earth's  surface  in  terms  of  the  lOP's  generated  ty  the 
system  [Eef.  1,  p.  14].  With  the  elimination  of  all  system¬ 
atic  cr  bias  errors,  the  terms  repeatability  and  predicta¬ 
bility  become  identical.  Hydrographic  surveyors  usually 
work  toward  this  condition,  although  it  is  not  always 
achievable. 

Heinzen  [Eef.  2]  and  Burt  [Hef.  3]  have  presented 
several  techniques  for  quantifying  the  repeatable  accuracy 
for  offshore  positions.  These  technigues  have  roots  in  the 
statistical  treatment  of  random  error.  Although  the  methods 
have  teen  well  documented,  no  single  criterion  to  classify 
the  accuracy  of  a  hydrographic  position  has  been  agreed  upon 
ty  the  international  hydrographic  community. 

Preceding  the  development  of  automation  in  hydrograpb 
data  acquisition  and  processing,  the  task  of  calculatirg  < 
accuracy  figure  to  attach  to  each  position  in  a  hydrograpti 
survey  was  unthinkable.  To  ensure  overall  accuracy  in  a 
survey,  certain  generalizations  were  developed  to  act  as 
guidelines.  For  example,  the  U.S.  Coast  and  Geodetic  Survey 
Hldrographic  Manual  [Eef.  4,  p.  217]  states  the  following 
concerning  the  strength  of  a  three-point  fix: 

The  fix  is  strong  when  the  sum  of  the  two  angles  is 
equal  to  or  greater  than  180°  and  neither  angle  is  less 
than  3C°.  The  nearer  the  angles  equal  each  other  the 
stronger  will  be  tie  fix. 

Generalizations  of  this  type  provided  useful  qualitative 
guidance  for  assuring  a  degree  cf  positional  accuracy  and 
many  are  still  in  existence  today. 

With  the  aid  of  computers,  the  hydrographer  now  has  the 
capacity  to  evaluate  the  accuracy  of  positioning  data  for  an 
entire  survey.  An  accuracy  figure  can  be  computed  for  each 
position  in  a  survey  and  stored  in  a  data  base  along  with 


ether  survey  information.  This  figure  may  provide  useful 
information  for  users  of  the  data,  as  well  as  a  yardsticK 
for  the  hydrographer  to  evaluate  the  quality  of  the  work. 
Furthermore,  a  presurvey  accuracy  analysis  enables  a  survey 
to  ce  designed  to  meet  desired  specifications. 


E.  ACCOEACY  STANDAEES  FOfi  HYDBCGRAPHIC  POSITIONING 


In  1S82,  the  International  Hydrographic  Organization 
(IHC)  putlished  new  recommendations  for  error  standards 
ccncernirg  the  accuracy  of  hydrographic  positions.  These 
standards  [Ref.  5]  are: 


The  position  of  soundings,  dangers  and  all  other  signif¬ 
icant  features  should  be  determined  with  an  accuracy 
such  that  any  probable  error,  measured  relative  tc  shore 
control,  shall  seldom  exceed  twice  the  minimum  plottarle 
error  at  the  scale  cf  the  survey  (normally  1.0  mm  on 
paper) .  It  is  most  desireable  that  whenever  positions 
are  determined  bv  the  intersection  of  lines  of  position, 
three  such  lines'be  used.  The  angle  between  any  pair 
should  net  be  less  than  30°. 


Most  statisticiars  define  the  term  "probable  error"  as 
that  error  occurring  at  the  50  percent  probability  level. 
However,  the  author  cf  the  IHC  standards.  Commodore  A . H . 
Cooper  RAN  (Ret.)  has  stated  that  the  term  "probable  error" 
was  interded  to  have  no  statistical  significance.  Munson 
interpreted  the  words  "shall  seldom  exceed"  to  mean  10 
percent  of  the  time  [Eef.  6],  Using  this  interpretation, 
the  first  sentence  of  the  specification  might  be  written: 

The  position  of  soundings,  dangers  and  all  other  signif¬ 
icant  features  should  be  determined  with  an  accuracy 
such  that  any  error  in  position  measured  relative  tc 
shore  control  will  fall  within  a  circle  with  radius  cf 
the  minimum  piottable  error  at  the  scale  of  the  survey 
(normally  1.0  mm.  on  paper),  with  90  percent  confidence. 

The  specification  in  this  form  could  be  evaluated  quantita¬ 
tively.  The  criterion  for  defining  accuracy  in  terms  cf  a 
fixed  probability  is  common  in  the  field  of  surveying.  For 
example,  the  standards  of  accuracy  developed  for  geodetic 


These  errors  are  usually  small  in  magnitude  and  can  le  elim¬ 
inated  by  proper  adjustment  of  the  instrument  by  either  the 
manufacturer  or  a  qualified  technician. 

The  field  hydrographer  has  ultimate  control  over 
the  geometric  systematic  errors  associated  with  a  theodo¬ 
lite.  In  range- azim u th  positioning  the  theodolite  and 
transmitter  may  occupy  the  same  horizontal  control  station. 
If  the  theodolite  is  not  set  directly  over  the  station  a 
resultant  systematic  error  will  occur  in  all  measurements. 

It  can  be  shown  that  these  errors  are  non-linear  but  do 
follow  a  mathematical  relationship.  Likewise,  if  the  trans¬ 
mitter  is  net  located  directly  over  the  station,  a  similar 
type  cf  bias  occurs.  Depending  on  the  eccentricity  cf  the 
theodolite,  the  vessel's  range  from  the  theodolite,  and  the 
scale  of  the  survey--these  errors  can  seriously  affect  the 
absolute  accuracy  of  the  offshore  positions. 

In  a  similar  fashion,  it  is  also  imperative  to 
position  the  target  directly  over  the  horizontal  control 
station  used  as  an  initial.  Failure  to  do  this  will  result 
in  an  error  which  will  be  propagated  to  offshore  positions. 

Many  situations  arise  in  the  field  where  it  is 
advantageous  to  set  a  transmitter  and  theodolite  over  a 
single  horizontal  control  station.  Frequently  it  is 
feasible  to  construct  a  platform  to  accommodate  both  instru¬ 
ments;  in  a  case  where  it  is  not,  the  position  of  an  eccen¬ 
tric  horizontal  control  station  near  the  original  station 
should  he  determined  and  that  station  used  for  the  location 
of  one  of  the  instruments.  The  theodolite  and  the  trans¬ 
mitter  then  occupy  two  known  stations  and  the  geometric 
source  of  systematic  error  is  eliminated. 

h.  Electronic  Ranging  Systems 

The  systematic  errors  associated  with  electronic 
positioning  systems  are  complex  in  nature  and  functions  of 


2.  Sys tema tic  Errors 

Systematic  errors  occur  with  the  same  sign,  usually 
of  similar  magnitude,  and  can  te  expressed  in  terms  of  a 
mathe iratical  model.  Systematic  errors  follow  a  defined 
pattern  and  occur  in  a  number  cf  consecutive  related  obser¬ 
vations.  Repetition  cf  measurements  does  nothing  tc  mini¬ 
mize  their  effect.  Ir  the  case  of  hydrographic  positioning, 
systematic  errors  are  identified  and  modeled  by  calibration 
of  the  measuring  instrument  against  a  known  standard.  Ihe 
following  is  a  brief  discussion  concerning  systeniatic  errors 
and  their  treatment  in  relation  to  hydrographic  positioning 
equip  me  r.t. 

a.  Theodolites 

In  nearshore  surveys  the  theodolite  is  used 
primarily  for  range-a2imuth  and  azimuth-azimuth  positioning. 
Systeniatic  errors  associated  with  the  theodolite  can  be 
classified  into  two  groups:  those  associated  with  the  phys¬ 
ical  design  of  the  instrument  and  those  involving  the  geom¬ 
etry  of  the  positioning  scheme.  Some  sources  of  systematic 
errors  [Ref.  8]  associated  with  the  physical  characteristics 
cf  a  theodolite  are: 

i.  The  horizontal  circle  may  be  eccentric, 

ii.  Graduations  on  the  horizontal  circle  may  not  be 
uniform. 

iii.  The  horizontal  axis  of  the  telescope  (about  which 
it  rotates)  may  not  be  perpendicular  to  the 
vertical  axis  of  the  instrument, 

iv.  The  longitudinal  axis  of  the  telescope  may  net  te 
normal  to  the  horizontal  axis, 

v.  The  telescope  axis  and  the  axis  of  the  leveling 
buttle  may  net  be  parallel. 
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range-azimuth  fix.  A  range  and  an  azimuth  are  generated 
from  a  known  control  station  to  the  vessel’s  position.  A 
second  control  station  is  used  to  fix  the  initial  azimuth;  a 
third  shore  control  station  is  located  TO  meters  from  the 
initial  station  and  its  coordinates  are  mistakenly  used  for 
the  initial  station  in  plotting.  The  resultant  hydrographic 
position  is  in  error,  but  this  error  will  not  be  easily 
distinguished. 

Although  most  blunders  have  their  origin  in  human 
carelessness,  some  can  be  attributed  to  equipment  malfunc¬ 
tion.  For  example,  microwave  systems  which  generate  IOF's 
are  known  to  become  unsteady  under  certain  conditions. 
Spuricus  range  readings  resulting  from  signal  reflections 
can  be  recorded  as  true  positioning  data.  Ir.  this  case,  the 
blunder  may  or  may  not  be  easily  detected. 

In  automated  data  acquisition  systems,  software  has 
teen  developed  tc  detect  the  occurrence  of  anomalous  range 
readings.  By  inputting  a  course  and  speed  of  a  vessel  trav¬ 
eling  along  a  line,  the  computer  can  determine  if  the 
recorded  position  is  valid  based  on  the  principle  of  dead 
reckoning.  If  the  recorded  position  is  found  to  be  invalid 
the  hydrographer  will  be  immediately  alerted  to  the  situ¬ 
ation  and  can  take  action  to  remedy  the  problem.  In  nor- 
automated  systems  the  principle  of  dead  reckoning  is  applied 
manually.  Given  the  course  and  speed  of  the  vessel,  the 
validity  of  the  position  can  be  checked  with  spacing 
dividers.  This  involves  checking  the  spacing  between  fixes 
recorded  before  and  after  the  position  in  question. 

Before  any  type  of  error  analysis  is  to  be  performed 
on  the  hydrographic  positioning  data,  it  is  essential  that 
all  blunders  be  identified  and  properly  treated.  In 
general,  careful  planning  coupled  with  thorough  checking 
will  micimize  the  occurrence  cf  blunders. 
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Figure  2.3  Geometry  of  an  Azimuth- Azimuth  Position 

Consider  the  following  as  an  example  of  a  blunder 
associated  with  range-range  geometry.  An  offshore  position 
is  to  be  determined  by  the  intersection  of  two  electronic 
LOP's  generated  from  transmitters  located  on  known  shore 
stations.  The  vessel  is  working  west  of  a  shoreline  that 
runs  generally  in  a  north-south  direction.  As  the  hydrc- 
grapher  faces  the  stations  from  sea,  the  southern  shore 
station  is  mistakenly  identified  as  left  and  the  northern 
shore  station  as  right.  The  resultant  offshore  position 
will  plot  to  the  east  of  the  base  line.  This  blunder  is 
readily  detected  and  can  be  easily  remedied. 

Sot  all  types  of  blunders  are  so  easily  detected. 
Suppose  an  offshore  position  is  to  be  determined  by  a 


22 


E.  CLASSES  OF  EHBOHS 

All  hydrographic  positioning  measurements  are  subject  to 
error.  Lhe  following  sections  discuss  categories  of  errors 
and  methods  used  to  treat  these  errors. 

1 .  Blunders 

Blunders  are  gross  mistakes  which  are  generally  due 
to  the  carelessness  of  the  observer.  Blunders  can  vary  in 
magnitude,  ranging  from  large  errors  which  are  easily 
detected,  to  small  errors  which  may  be  barely  distinguished. 
They  can  be  detected  by  making  repeated  observations  cr  by 
carefully  checking  the  data  in  the  processing  phase. 

Blunders  occur  in  various  forms  and  most  can  be  avoided  by 
carefully  planning  tte  data  acquisition  process. 
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in  this  arrangement  tut  systems  employing  a  laser  can  also 
be  used  for  short-range  work.  Another  LOP  is  generated  by 
fixing  an  azimuth  from  a  shore  control  station  to  the 
vessel.  A  second  control  station  is  used  for  an  initial 
azimuth  by  the  observer.  Azimuth  determinations  can  be  made 
after  observing  directions  with  a  theodolite  as  an  observer 
tracks  the  moving  vessel. 

There  are  t wc  ways  to  determine  a  range-azimuth 
position.  The  most  common  way  is  to  have  the  theodolite  and 
the  transmitter  occupy  the  same  shore  control  station. 

Hence,  the  angle  of  intersection,  g,  of  the  LOP's  is  always 
50°.  This  arrangement  is  commonly  used  by  the  National 
Ccean  Service  (NOS)  for  large-scale  nearshore  surveys. 

The  other  way  is  to  have  the  theodolite  and  the 
transmitter  occupy  two  different  control  points.  Then  the 
geometry  is  similar  to  that  of  the  range-range  position. 

The  angle  of  intersection,  g,  is  computed  by  trigoncme t ric 
relationships  among  the  azimuth  of  a  line  between  the  shore 
stations,  the  observed  direction  to  the  vessel,  and  the 
measured  range  to  the  vessel. 

4 .  Azimuth- Azimuth 

Azimuth-azimuth  positioning  geometry  is  used  for 
nearshore  high-accuracy  surveying.  Theodolites  are  set  over 
two  control  stations  on  shore.  The  vessel  is  sighted  on 
simultaneously  by  the  two  theodolite  observers,  generating 
two  visual  LOP’s  whose  intersection  define  the  vessel's 
location.  Initial  azimuths  are  fixed  by  sighting  on  control 
stations  which  are  visible  to  the  observers. 

The  angle  of  intersection  for  an  azimuth-azimuth 
position  is  dependent  on  the  geometric  relationships  between 
the  occupied  stations,  the  initial  stations,  and  vessel's 
position  (Fig.  2.3).  Assuming  that  theodolite  observers 
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where  the  term  1/sin  (a/2)  is  called  the  lane  expansion 
factor.  The  angle  of  intersection  ,B ,  between  the  twc 
hyperbolas  is  then  given  by 


Figure  2.2  Geometry  of  a  Hyperbolic-Hyperbolic  Position 

3  .  Range- A zia at  h 

This  positioning  geometry  is  used  for  nearshore, 
line-cf-sight  surveys.  One  LOP  is  generated  by  an  elec¬ 
tronic  range  originating  from  a  transmitter  located  on  a 
shore  control  station.  A  microwave  system  is  commonly  used 


Hyperbolic  location  methods  can  be  divided  into  two 
groups  based  on  the  electronic  principles  used  to  define  the 
distance  differences  £Ref.  7,  p.  87],  Loran  is  an  example 
of  a  pulse  system  in  which  the  differences  in  times  of 
arrival  of  pulses  transmitted  ty  the  master-slave  combina¬ 
tions  are  translated  into  distance  differences.  The  resul¬ 
tant  position  has  no  lane  ambiguity  and  is  easily  resolved. 
The  second  method  of  hyperbolic  positioning  involves  meas¬ 
uring  a  phase  differerce  from  two  master-slave  combinations 
at  the  vessel's  position.  The  phase  difference  translates 
into  a  fractional  lane  count  which  in  itself  provides  an 
ambiguous  position.  This  ambiguity  is  resolved  ty  using  a 
whole-lane  counter  which  is  initialized  at  a  known  geograph¬ 
ical  point.  In  hyperbolic  positioning,  the  ship  is  in  a 
passive  mode  and  the  system  can  be  used  by  many  vessels. 

The  angle  of  intersection  between  the  two  hyperbolas 
can  be  computed  by  first  defining  the  following  quantities: 
is  the  length  of  red  base  line, 

S  is  the  length  of  green  base  line, 

R  is  the  distance  between  master  and  vessel's  position  ?, 
m  c 

Rr  is  the  distance  from  red  slave  to  point  P, 

R  is  the  distance  from  green  slave  to  point  P, 

s 

is  the  angle  between  lines  P!1  and  PR,  and 

a  is  the  angle  between  lines  PM  and  PG. 

g 

The  spacing  between  lanes  increases  with  distance 
from  the  master-slave  pair.  The  lane  widths  along  the  base 
line  are 


w'  A 
r  =  r 


and 


w'  A 
9  *  J. 

r 


(2.5) 


Then  the  lane  widths  at  any  point  P  are 

A  -  X. 


w  =  r  ,  1  a  and  w  =  g  ,  _ l  \ 

r  r  (  sin  [aj2)  ]  9  2  (  sin  (a  / 2) 


(2.6) 
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Figure  2.1  Geometry  of  a  Range-Range  Position 

2 .  Hyperbolic-H  vcerbolic 

Hydrographic  positioning  by  hyperbolic-hyperbolic 
geometry  utilizes  the  intersection  of  two  hyperbolas  each 
generated  about  a  pair  of  shore  control  stations.  A  hyper¬ 
bola  is  the  locus  of  points  in  which  the  difference  cf 
distance  from  two  fixed  points  is  always  constant.  A  three- 
station  hyperbolic  net  is  the  most  commonly  used  hypertclic 
mode  for  offshore  survey  (Fig.  2.2).  One  family  of  hyper¬ 
bolas  (Red)  are  generated  about  a  master  station,  M ,  and  a 
slave,  R;  while  a  second  family  of  hyperbolas  (Green)  are 
generated  with  respect  to  the  master  and  a  second  slave,  G. 
For  the  first  family  cf  hyperbolas,  the  control  points  K  and 
R  act  as  the  foci,  while  points  M  and  G  act  as  the  foci  for 
the  second  family. 


problem  cf  lane  ambiguity  must  be  addressed.  Ranges  are 
expressed  in  full  and  partial  lane  counts  where  a  lane  width 


w  =  ! 
w  2 


(2.2) 


where  A  is  the  wavelength  of  the  transmitting  frequency,  f, 
and  given  by 


(2.3) 


Medium- range  systems  commonly  in  use  today  are  Cubic 
Western's  "ARGO,"  Hasting  Raydist's  "Raydist,"  and  Odom 
Cffshcre's  "Hydrotrack. '• 

The  angle  of  intersection  associated  with  a  range- 
range  position  is  computed  frcm  a  simple  trigonometric  rela¬ 
tionship.  The  vessel's  position  P  (Fig.  2.1)  is  determined 
by  the  intersection  cf  the  ranges  from  the  left  and  right 
shore  stations,  R1  and  R2  respectively.  3  is  the  base  line 
distance  computed  between  the  two  known  shore  stations. 

Since  the  range  circles  from  the  shore  stations  intersect  at 
two  points,  it  is  necessary  for  the  plotter  to  recognize 
which  side  of  the  base  line  the  vessel  is  on  in  order  to 
eliminate  the  ambiguity.  The  angle  of  intersection  cf  the 
two  LOP's  (8)  is  given  by  the  law  of  cosines 


8  »  180°  -  Arc  cos  (  ) 


(2.4) 


In  qualitative  terms,  the  fix  is  strongest  when  6  approaches 
900.  Most  hydrographic  specif ications  limit  the  angle  cf 
intersection  from  a  ninimum  of  30°  to  a  maximum  of  150°. 
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An  electronic  positioning  system  may  be  active  cr 
passive.  In  an  active  system,  a  transmitter  from  the  survey 
launch  keys  the  transmission  of  ranges  from  the  shore 
station.  In  turn,  the  signals  generated  from  the  shore 
stations  (slaves)  are  then  received  by  the  launch.  An 
active  system  is  limited  to  a  finite  number  of  users, 
usually  not  more  than  about  four-  The  number  of  users  cf  a 
passive  system  is  unlimited  as  the  survey  launch  requires 
only  a  receiver  which  is  constantly  listening  for  signals 
which  are  being  transmitted  from  shore. 

Short-range,  cr  line-of-sight,  positioning  systems 
are  used  for  nearshore  hydrographic  surveys.  These  systems 
operate  in  the  microwave  region  of  the  electromagnetic  spec¬ 
trum  (3  to  10  GHz).  A  distance  is  determined  by  observing 
the  time  needed  for  a  pulse  to  travel  from  a  master  tran¬ 
sponder  located  aboard  the  survey  vessel  to  a  remote  tran¬ 
sponder  on  shore  and  back  to  the  master  transponder. 

Knowing  the  average  velocity  of  the  electromagnetic  pulse, 
the  distance  D  is  then 


where  c  is  the  group  velocity  of  the  wave  packet  and  t  is 
the  two-way  travel  time.  Short-range  systems  which  are  in 
wide  use  today  are  Eacal  Decca's  "Tr isponder"  and  Motorola’s 
"Mini-Banger.’’  These  systems  have  direct  range  readout  and 
are  readily  interfaced  into  a  navigational  computer  and  a 
data  acguisition  system.  Both  systems  are  active  and  user 
limited. 

Medium-range  positioning  systems  operate  in  the  1- 
to  5—  KHz  frequency  range  of  the  electromagnetic  spectrum.  A 
distance  is  determined  by  measuring  the  phase  relationship 
between  transmitted  and  received  waves.  These  systems  are 
usually  referred  tc  as  continuous  wave  systems  and  the 


II.  8AT0RE  OF  THE  PROBLEM 

The  development  cf  an  accuracy  figure  for  offshore  posi¬ 
tions  is  inherently  tied  to  the  geometry  of  the  positioning 
method  and  the  errors  which  are  associated  with  the  posi¬ 
tioning  eguipment  that  is  used.  This  chapter  will  discuss 
the  geometric  and  statistical  elements  involved  in  deter¬ 
mining  an  offshore  position  and  presents  several  methods  for 
guantifying  repeatable  accuracy. 

A.  HIDECGRAPHIC  POSITIONING  GEOMETRIES 

An  offshore  fix  can  be  determined  by  the  intersection  cf 
two  or  mere  LOP's.  These  LOP's  may  be  generated  by  elec¬ 
tronic  or  visual  means.  Working  toward  the  development  of 
an  accuracy  index,  it  will  be  necessary  to  compute  the  ancle 
of  intersection  of  the  LOP's  associated  with  different  posi¬ 
tioning  geometries.  The  following  sections  discuss  the 
geometry  cf  conventional  offshore  positioning  methods  and 
ways  to  compute  the  angles  of  intersection.  This  thesis 
will  not  address  the  geometry  involved  in  a  three-point 
sextant  fix. 

1 .  Range-Range 

Establishing  an  offshore  fix  by  range-range  geometry 
involves  measuring  distances  electronically  from  fixed  posi¬ 
tions  on  shore  to  the  vessel's  unknown  location.  Ranges  can 
be  determined  by  measuring  the  elapsed  time  between  trans¬ 
mission  and  receipt  of  a  radio  pulse  or  by  comparing  the 
phase  of  the  transmitted  wave  with  the  phase  of  the  received 
wave  [Ref.  2].  In  each  case,  transmitters  are  set  on 
stations  on  shore  whose  coordinates  are  determined  by 
precise  land  survey  methods. 
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method  of  classification  is  a  useful  index  for  quantifying 
the  accuracy  of  positions.  The  computed  radii  of  the  9C 
percent  confidence  circles  can  serve  as  an  accuracy  figure 
that  can  be  attached  to  each  position  in  a  survey  and  stored 
in  a  cata  base. 

The  third  objective  of  this  thesis  is  to  demonstrate 
that  a  presurvey  analysis  can  be  used  in  designing  posi¬ 
tional  accuracy  to  meet  specifications.  The  existing 
general  guidelines  for  planning  can  be  tetter  defined.  For 
example,  in  planning  a  survey  hydrographers  usually  lay  out 
circles  which  delimit  the  30°  and  150°  boundaries  that 
define  the  minimum  and  maximum  allowable  intersection  angles 
between  two  LOP’s.  As  a  means  to  meet  accuracy  require¬ 
ments,  it  can  be  shown  that  these  limits  should  vary  based 
on  the  scale  of  the  survey  and  the  precision  of  the  posi¬ 
tioning  equipment. 
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control  surveys  have  their  origin  in  probability  theory. 
Procedures  for  obtaining  first-order  geodetic  positions 
require  sixteen  repeated  theodolite  observations  of  each 
direction.  Lower  order  positions  require  fewer  numbers  of 
observations.  Given  the  precision  of  one  observation  of 
each  direction,  it  car  be  demonstrated  that  increasing  the 
number  of  observations  coincides  with  increasing  the  prob¬ 
ability  of  the  direction  falling  within  specified  limits. 

Eegarding  accuracy  determinations,  there  are  several 
problems  unique  to  hydrographic  surveying.  Whereas  stan¬ 
dards  for  other  types  of  surveys  rely  on  multiple  observa¬ 
tions  of  the  same  quantity,  the  accuracy  of  a  hydrographic 
position  must  be  evaluated  in  terms  of  a  single  observation 
(which  may  be  the  intersection  of  two  or  more  LOP’s). 

Diverse  methods  for  obtaining  a  hydrographic  position  exist 
and  these  methods  must  all  be  evaluated  using  the  Same 
criterion.  Also,  there  is  a  broad  spectrum  of  equipment 
used  in  hydrographic  positioning  and  in  many  cases  the 
precision  of  this  equipment  is  not  well  defined. 

C.  CEJECTIVES 

A  need  exists  to  give  quantitative  meaning  to  the  accu¬ 
racy  specifications  set  forth  by  the  IHO.  One  of  the  objec¬ 
tives  of  this  thesis  is  to  demonstrate  that  defining  the 
specifications  in  terms  of  the  fixed  90  percent  confidence 
level  is  a  valid  interpretation.  By  defining  what  the  spec¬ 
ifications  imply,  procedures  can  be  developed  to  meet  the 
stand  ards. 

A  second  objective  of  this  thesis  is  to  apply  the  theory 
of  errcrs,  associated  with  hydrographic  positioning,  to  a 
data  set.  This  analysis  involves  classifying  positioning 
data  acquired  in  a  survey  based  on  the  radii  of  circles  of 
equivalent  probability.  It  will  be  demonstrated  that  this 
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many  variables.  Munson  [Ref.  9,  p.  4]  addresses  several 
problems  associated  with  short-range  systems  used  i  r.  hydro¬ 
graphic  surveys.  The  most  common  problems  with  short-range 
systems  are  variation  in  range  and  calibration  drift  with 
time.  Variations  in  internal  eguipment  time  delays  in  the 
transmitter,  the  transponder,  or  the  receiver  can  induce 
errors  in  measured  ranges.  For  pulse  systems  such  varia¬ 
tions  can  occur  due  to  temperature  dependence  of  components 
and  fluctuations  in  signal  strength  at  the  transponder. 
Multipath  effects  are  also  a  problem.  Under  some  circum¬ 
stances  a  reflected  wave  and  the  directly  transmitted  wave 
arrive  with  a  phase  difference  of  180°.  Cancellation  or 
fading  of  the  directly  transmitted  signal  can  result. 

NOS  conducts  base  line  calibrations  of  short- 
range  positioning  systems  periodically  during  the  course  of 
a  survey  to  minimize  or  eliminate  systematic  error.  In  this 
process,  a  transmitter  and  receiver  are  each  placed  ever 
control  stations  on  shore  and  the  measured  range  is  compared 
to  the  true  range.  In  this  way  the  systematic  error  is 
eliminated  by  zeroing  the  instrument  or  by  applying  a 
constant  correction  to  raw  data.  System  checks  are 
performed  daily  to  assure  there  is  no  drift  from  the  orig¬ 
inal  calibration.  A  check  can  be  accomplished  by  comparing 
a  positicn  defined  by  the  ranging  system  to  a  known  fixed- 
point  position,  to  a  sextant  fix  position,  or  an  intersec¬ 
tion  position. 

Munson  [ Bef .  9,  p.  5]  also  discusses  sources  of 
systematic  errors  associated  with  medium-range  systems.  The 
most  significant  systematic  errors  occur  as  a  function  of 
position  due  to  varying  propagation  velocity.  The  medium- 
range  electronic  signal  propagation  velocity  depends  on  the 
surface  conductivity  and  transmission  path  (over  water,  over 
land,  or  over  different  types  of  land) .  3ecause  of  this 
dependence,  systematic  errors  as  a  function  of  position 


occur  at  different  effective  phase  velocities.  Knowing  the 
propagation  velocity  to  use,  or  the  phase  correction  tc  make 
as  a  function  of  range,  is  a  problem.  Sky  wave  and  storm 
interference  also  pose  problems.  At  extreme  ranges  of  oper¬ 
ation,  sky  wave  interference  can  affect  the  more  predictable 
ground  wave,  especially  during  nighttime  operations.  Lane 
ambiguities  are  also  a  problem.  Most  systems  are  inherently 
ambiguous  and  must  be  zero  set  and  continually  monitored  for 
lane  jumps  or  loss  of  signal  which  results  in  the  loss  of 
lane  count. 

NOS  uses  several  techniques  to  determine  the 
systematic  error  associated  with  medium-range  positioning 
systems.  These  techniques  involve  determining  a  whole  and 
partial  lane  count  for  phase  comparison  systems.  Two  of  the 
more  widely  used  techniques  are  comparison  of  three-point 
sextant  fix  positions  to  positions  determined  by  the  elec¬ 
tronic  ranging  system  and  calibration  of  the  electronic 
system  at  a  fixed  point.  In  both  techniques  the  whole  lace 
counts  are  fixed  by  the  calibration;  correctors  to  the 
partial  lace  count  are  determined  and  applied  to  the  raw 
ranging  data. 

3  •  Fandom  Errors 

Fandom  errors  are  chance  errors,  unpredictable  in 
magnitude  or  sign,  and  are  governed  by  the  laws  of  prob¬ 
ability  [Ref.  10,  p.  1206],  They  are  errors  which  remain 
after  blunders  and  systematic  errors  have  been  removed. 
Random  errors  result  from  accidental  and  unknown  combina¬ 
tions  of  causes  and  ate  beyond  the  control  of  the  observer. 
Greenwalt  [Ref.  12,  p.  2]  states  they  are  characterized  by: 

i.  Variation  in  sign;  positive  errors  occur  with 
equal  frequency  as  negative  ones, 

ii.  Email  errors  occur  mere  frequently  than  large  errors 

iii.  Extremely  larce  errors  rarelv  occur. 


Random  errors  are  unique  to  specific  types  of  posi¬ 
tioning  equipment  and  vary  in  magnitude  depending  on  the 
precision  of  the  instruments  that  are  used.  The  following 
section  outlines  statistical  methods  for  their  treatment. 

C.  TSI ATHENT  OF  RAN IC9  ERRORS. 

1 .  Cne-Di mensional  Errors 

Certain  basic  statistical  quantities  must  first  be 
defined  in  the  analysis  of  random  errors.  Consider  a  vessel 
moored  securely  to  a  fixed  offshore  platform.  A  number  of 
ranges,  n,  from  a  microwave  transmitter  located  on  a  shore 
control  station  are  recorded.  The  mean  of  these  observa¬ 
tions  is 

n  x, 

U  =»  E  — 

x  i«l  n  (2.  1 1) 

where  x^  represents  an  individual  observation.  The  standard 
error,  s,  of  the  observations  is  then 

(2.  12) 

where  the  quantity  (x^  )  is  referred  to  as  the  residual, 

cr  true  error,  v^ ,  of  a  particular  observation.  As  n  gets 
very  large,  the  factor  1/n  can  be  substituted  for  1/(n-1)  in 
Equation  2.12.  likewise,  in  treating  the  large  sample,  o 
can  be  substituted  for  s  and  m  for  ,  where  y  and  a  are 
the  mean  and  standard  error  of  the  entire  population. 

It  is  of  interest  to  determine  the  probability  of 
occurrence  of  a  particular  observation.  The  normal  cr 
Gaussian  distribution  equation  relates  the  residual  of  a 
particular  random  variable  with  the  probability  of  its 
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occurrence,  and  is  given  by 

P(V)  » 


(2.  13) 


a  Vr2ir 

The  plot  of  this  equation  yields  the  normal  distribution 

curve  (Fig.  2.4).  The  height  cf  the  curve  above  the 

vertical  axis  is  proportional  to  the  probability  of  a 

particular  error  occurring. 

The  probability  of  a  residual  falling  between  any 

two  residuals  v  and  v  can  be  computed  by  integrating 
12 

Equation  2.13  as 


P(v) 
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dv 


(2. 14) 


Figure  2.4  The  normal  Distribution 

This  integral  is  difficult  to  evaluate  analytically 

so  tables  have  been  compiled  to  aid  in  computations.  Fcr 

v  =  *a  and  v  =  -a,  it  can  be  shown  that  P(v)  =  0.6827.  In 
1  2 

ether  words,  the  probability  that  a  particular  observation 
will  fall  within  ♦  1o  of  the  mean  is  68.27  percent. 


Fetarning  to  the  example  of  the  vessel  moored  tc  the 
cffshcre  platform,  the  mean  and  the  standard  error  for  the 
observations  are  easily  computed.  Kith  this  information  and 
Equation  2.14,  the  probability  of  a  range  error  falling 
within  specified  limits  can  be  computed.  Conversely,  by 
fixing  a  probability,  the  associated  limits  of  the  range 
error  can  be  computed  .  In  statistical  terms,  a  particular 
observation  will  fall  within  specified  limits  with  a  certain 
confidence. 

Actual  values  of  one-dimensional  standard  errors  for 
hydrographic  positioning  equipment  are  a  subject  of  debate 
betweec  manufacturers  and  users.  Some  manufacturers  of 
microwave  positioning  equipment  claim  standard  errors  of  +1 
meter.  Cn  the  other  hand,  Munson  [Ref.  9,  p.  6]  states  that 
microwave  systems  demonstrate  accuracies  of  3  meters  at 
short  ranges  but  show  larger  errors  at  ranges  of  15  km  and 
greater.  NOS  assumes  a  3-meter  standard  error  in  all  of  its 
short-range  accuracy  computations.  It  is  apparent  that 
further  study  is  needed  to  adequately  define  the  nature  of 
errors  associated  with  electronic  positioning  equipment. 

Kaltz  [Ref.  13]  performed  an  extensive  study  to 
determine  the  pointing  error  of  a  Wild  T-2  theodolite.  His 
results  showed  that  the  pointing  error  associated  with  this 
instrument  under  hydrographic  survey  conditions  was  about 
1.3  meters  and  was  independent  of  distance. 

2 •  Iwo-Di mensional  Errors 

Ihe  intent  of  this  paper  is  to  apply  statistical 
methods  developed  by  ethers  tc  a  hydrographic  data  set 
containing  two-dimensional  errors  which  are  defined  by  two 
random  variables.  Lengthly  and  complex  derivations  are  not 
presented.  Burt  [Bef.  3]  and  Heinzen  [Ref.  2]  show  adequate 
derivations  of  formulas  associated  with  two-dimensional 
errors  ard  can  be  referenced  for  full  details. 
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The  following  assumptions  are  made  concerning  two- 
dimensional  errors  associated  with  intersecting  LCP's: 

i.  The  random  enters  of  each  LOP  are  normally 
distributed. 

ii.  Systematic  or  bias  errors  have  been  removed  from 
the  observations. 

iii.  The  intersecting  LOP's  are  coplanar. 

iv.  The  error  LOP’s  are  parallel  to  the  exact  LOP's. 


In  developing  a  usable  mathematical  model  for  accuracy 
determinations,  the  four  assumptions  hold  to  a  high  degree 
for  all  hydrographic  positioning  geometries. 

Consider  again  the  vessel  moored  to  a  fixed  offshore 
platform.  Assume  two  ranges  are  measured  from  two  different 
shore  control  stations  at  the  same  time  and  that  the  range 
readings  are  uncorrelated.  The  observation  of  this  pair  of 
ranges  is  repeated  many  times.  After  a  large  number  of 
observations,  the  means  and  standard  errors  of  the  indi¬ 
vidual  ranges  are  determined.  Suppose  the  mean  ranges,  or 
the  actual  LOP's,  intersect  at  an  angle  of  90°  and  that  the 

computed  standard  errors  are  egual  (a  -a  }-  if  each  data 

1  2 

pair  (x^ ,y  )  is  plotted,  the  spread  of  points  about  the  mean 
coordinates  results  in  a  circular  cluster  (Fig.  2.5).  A 
higher  density  of  points  occurs  near  the  intersection  of  the 
mean  ranges  and  the  density  of  points  decreases  outward  from 
the  intersection  of  the  mean  ranges. 

In  this  special  case,  which  is  called  a  circular 
normal  distribution,  the  probability  of  a  point  falling 
within  a  specified  radius,  R,  from  the  intersection  of  the 
mean  ranges  is 


P(R)  ■  1  -  e 


l  Rai 


3  ^ 


(2.  15) 


where  o  =  a  =  a_  and  is  defined  as  the  circular  standard 
1  2  C 

error.  Using  Eguaticr  2.15,  R  can  be  computed  by  fixing 

P  (R)  ,  cr  ccnversly,  P  (R)  can  be  computed  by  fixing  R. 

Letting  R  =  o  -a  -a  ,  then  P(R)  =  0.  3935.  In  other 
1  2  c 

words,  39.35  percent  cf  all  errors  in  a  circular  normal 
distribution  are  not  expected  to  exceed  the  circular  stan¬ 
dard  error  [Ref.  12,  pp.  25-26  ]. 


Figure  2.5  Circular  Normal  Distribution 


In  the  case  where  the  two  uncorrelated  LCP’s  inter¬ 
sect  at  an  angle  other  than  90°  or  cr  ? a  ,  the  contours  cf 

l  2 


equal  density  are  ellipses  centered  about  the  point  defined 
by  the  intersecting  ICP's  (Fig.  2.6) .  The  two-dimensional 
probability  density  function  becomes  [Ref.  1,  p.  136] 


P(v  ,v  ) 
'  x*  y 


_ 1_ 

2irc  o" 


K2 


(2.  16) 


x  y 


Figure  2.6  Error  Ellipse  Formed  by  Two  Dncorrelated  LCP's 
where 

7  is  the  residual  in  the  direction  of  the  semi-major  axis 

x 

of  the  error  ellipse. 


7  is  the  residual  in  the  direction  of  the  semi- mine r 


ax  is  the  standard  error  in  the  direction  of  the  semi¬ 
major  axis, 

a  is  the  standard  error  in  the  direction  of  the  semi- 

y 

miner  axis. 


v2  v2 


k*  =  -4  +  4 

N  a 1  a 
x  y 


(2.  17) 


The  solution  of  Eguation  2.16  with  values  of  K  for  different 
P's  yields  the  results  in  Table  I  [Bef.  12,  p.  23].  For  a 
39.35  percent  probability,  the  axes  of  the  ellipse  are 
1.0000  a  and  I.OOoO  a  ;  for  a  50  percent  probability,  the 


axes  are  1.1774 


id7  1.  1774  a  - 

y 


TABLE 

I 

Values  of  the 

Constant  K 

PROEAEILITY 

K 

39-35% 

1.0000 

50.00% 

1  .  1774 

63.21% 

1. 4  1  42 

90.00% 

2.1460 

99.00% 

3.0349 

99.78% 

3.5000 

The  error  ellipse  can  be  used  for  accuracy  computa¬ 
tions  by  developing  relationships  for  o  and  o  in  terms  cf 

x  y 

the  initial  information  a  ,  a  ,  and  6  *  3owditch  [Bef.  10, 

1  z 

p.  1213]  gives  the  following  equations  for  independent  top's 
relating  these  quantities: 


2sin'S 


{  o2  ♦  o2  + 
1  2 


y  +  o\ 

l  2 


fsTnr6-oTo7’  } 
1  2 


(2. 18) 
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and 


,2  = 


I 


1 


sin  3 


{  a2  +  a*  -  V  lo*  +  a2)*  -  4sin‘3  aioi 
12  12  12 


[2.  19) 


In  these  equations,  g  is  assumed  to  be  the  acute  angle 
between  the  LOP's. 

In  certain  special  cases,  the  above  equations  take 
on  more  manageable  fcrms.  In  range-range  and  azimuth- 
azimuth  positioning  it  is  often  assumed  that  a ^  ~  a  =  o. 

Equations  2-18  and  2.19  then  reduce  to 

fZ  _ 

°x  *  2s in ( ^8 )  ^°) 

and 

-  ^  - 

ay  2cos(48)  (2-21) 


In  the  concentric  range-azxmuth  case,  c  *  a  f  and  B 

1  2 

equals  90°.  Equations  2.18  and  2.19  then  simplify  to 


x  = 


(2.  22) 


and 


o 


y 


o 

2 


(2.23) 


where  a  >  a  and  a  >  a  • 

12  x  y 

The  case  for  correlated  LOP's  is  more  complex.  The 

calculation  of  o  and  o  involves  a  coordinate  transforms- 

x  y 

tion  from  a  linear  skewed  coordinate  system  to  an  unccrre- 
lated  rectangular  coordinate  system.  The  following 
discussion  is  taken  from  Heinzen  [Ref.  2,  pp.  49-53], 

Assume  a  hydrographic  position  is  established  by  the 
intersection  of  two  correlated  LOP's  (Fig.  2.7a).  LOP  1  and 
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*.  -  - - - - 1 

F .  ;  ,  i  *  '  ::  in:i.‘cr2dtions  for  Correlated  10  Pfs 

-  1  .  s :  »  •  *  ;  i*>  u*  .■  n.  the  shewed  coordinate 

l".  *. «  i ,  «.•  •  i  •  . :  its  -  and  o  .  The  semi-major  and 

1  2 

t.~a  :  id'.,  .i  t  . »  t-i  .  u  »  1  lip se  dre  not  coincident  with 

t  :.e  ikt«>  :  coo:  ic.jt.  s  •  -  tea  The  correlation  coeffi¬ 

cient  b  t  f  w e  c  :i  *■’. e  t  j  c  ^ <.  P  9  s  in  p  •  A  s  s  u  m  e  rj  ^  q  . 

‘  2  12 

The  st  tr.oard  errors  arc  correlation  coefficient  in  a 
correlated  rectangular  coordinate  system  with  axes  A  and  B 
must  new  be  determined.  A  coordinate  transformation  from 
the  skewed  system  to  the  correlated  rectangular  system  must 
be  made  yielding  the  standard  errors  along  the  new  coordi¬ 
nate  axes  (Fig.  2.7b) 

+  2d  a  a  cos8  +  a*)  "  °2 
a  Sin  B  1  12  1  2  2  2 
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and 


ab  3  a2 


(2.25) 


Ihe  correlation  coefficient  in  the  correlated  rectangular 
system  is 

)  {  1  +  P  tr^cosS  +  (^-)2 

1  2 


P  .  *  (  r2-  cosS  +  o  )  {  1  +  0  (r^)cosB  +  (-^cos^}"1*  _ , 

ab  o  12  12  o  o  2.zo) 

i  li 


To  determine  a  and  0  ,  a  second  coordinate  transformation 
x  y 

aust  be  performed  frca  the  correlated  rectangular  system  to 
an  uncorrelated  rectangular  system  with  axes  X  and  Y  (rig. 
2.7c).  The  semi-majcr  and  semi-minor  axes  of  the  error 
ellipse  are  then 


i  ♦ 


11  -  Pjb) 

( «i  ♦Tir 


(2.27) 


an  1 


+  <  - 


( 2 .  28 ) 


When  °12  =  0,  these  ecuations  become  identical  to  the 
simplified  versions  in  Bowditch  £Ref.  10]. 

The  orientation  of  the  semi-major  and  semi-miner 
axes  relative  to  the  intersecting  LOP's  is  the  third  param¬ 
eter  which  fixes  the  error  ellipse.  The  angle  9  (Figs.  2.6 
and  2.7)  is  measured  counter-clockwise  from  LOP  1  to  the 
semi-aajcr  axis  of  the  error  ellipse  £Bef.  11]  and  is  given 
by 

a2sin(2B)  +  2o  a  a  sln(S) 


1 


12  1  2 


9  »  j  arctan  {  — -  - 

c  o2cos(2B)  +  2o  ao  cos(S)  +  0 


} 


(2.29) 


12  1  2 


For  the  special  case  cf  a  -  a  and 

1  2 


P  =  0 , 
1  2 


••I 


(2.  20) 
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Ihe  orientation  of  the  error  ellij.se  in  an  orthogonal  coor¬ 
dinate  system  can  be  represented  by  adding  or  subtracting  9 
to  the  orientation  of  LOP  1.  Care  must  be  taken  on  deter¬ 
mining  the  quadrant  of  the  outcome-  As  a  general  rule,  the 
error  ellipse  always  lies  within  tae  acute  angles  formed  by 
the  intersecting  LOP's. 

Ihe  orientation  and  dimensions  of  the  error  ellipse 

provide  a  useful  index  for  evaluating  the  accuracy  of  a 

hydrographic  position.  Its  greatest  attribute  is  that  it 

accurately  represents  the  error  distribution  about  the 

intersection  of  two  ICP's  in  terms  of  a  fixed  probability. 

It  is  interesting  to  examine  the  variation  in  the  relative 

dimensions  and  orientations  of  error  ellipses  as  they  vary 

in  a  range-range  configuration  with  a  =  o  =  a  (Fig*  2.8)  . 

1  2 

Ihe  dimensions  of  the  ellipses  are  specified  by  Eguaticrs 

2.20  and  2.21  and  a  and  a  are  functions  of  g  only  for 

x  y 

fixed  a  .  Therefore,  the  dimensions  or  the  ellipses  remain 
constant  along  a  contour  of  constant  g  ;  only  the 
orientation  changes.  A  line  cf  constant  8  is  a  circle 
which  includes  stations  L  and  B.  Note  that  the  dimensions 
of  the  ellipses  for  8  's  of  30°  and  150°  are  identical.  The 
ellipses  about  the  9C°  angle  of  intersection  contour  are 
circles  and  represent  the  strongest  possible  positions  in 
this  scheme.  With  varying  8 ’s,  the  directional  nature  cf 
the  distribution  can  be  noted. 

3.  Circular  Precision  Indexes 

Although  the  error  ellipse  gives  a  true  representa¬ 
tion  cf  the  error  distribution  about  a  hydrographic  posi¬ 
tion,  its  use  has  certain  drawbacks.  The  cbaracterist ics  of 

the  ellipse  must  be  s  ecified  by  the  tnree  quantities  o  , 

x 

1  ,  and  9.  A  single  figure  for  evaluating  the  positional 

accuracy  cannot  be  used.  Greenwalt  £Ref.  12,  p.  26]  states 

that  when  cj  and  0  are  not  equal,  a  circular  error 

x  y 
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TAE1E  V 

Coordinates  of  Control  Stations 


S1ATICN  MilS 
USE  MON 

MUSSEL 

EEACH  LA3 


GJCEETIC  CCCRD. 

36°  36*  04.685"  N 
121°  52'  35.900"  W 

360  37*  18.151"  N 

1  2  10  54’  1  1.628"  Vi 

36°  36'  05.57  1"  A 
1210  52'  33.427"  W 


Mid  COORD  . 

Y  =  1992.  43  Di . 
X  =  4853.36  I . 

Y  =  4247.  42  m. 
X  =  2474.75  id. 

Y  =  2009. 86  m. 
X  =  4914.75  ni . 


B.  ACCURACY  ANALYSIS  OF  HYDROGRAPHIC  POSITIONING  DATA 

lie  objective  of  this  section  is  to  illustrate  how  the 
accuracy  of  hydr ogra £ hie  positioning  data  can  be  classified 
using  Eurt's  method  of  circles  of  equivalent  probability. 
The  radius  of  the  90  percent  confidence  circle  was  computed 
for  each  position;  it  provides  a  quantitative  measure  of 
repeatable  accuracy. 

Fcr  subsequent  accuracy  computations,  the  following 
assumptions  were  made: 

i.  The  standard  error  for  the  microwave  ranging 
system  used  in  the  range-range  and  range- 
azimuth  computations  is  3  meters, 
ii  .  For  azimuth-azimuth  and  range-azimuth  positions, 
the  pointing  error  of  the  theodolite  is  1.3 
meters  at  all  ranges. 

iii.  The  two  LOP's  involved  in  all  types 

cf  positioning  are  independent  (  p  =0). 

1  2 

iv.  The  data  are  tree  of  systematic  errors. 

Raw  range  and  azimuth  data  were  hand  logged  into  a  data 
file  for  processing.  A  modification  of  program  UCOMPS  was 
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IV.  RESULTS  AND  DATA  ANALYSIS 

A.  DATA  PROCESSING 

Automated  processing  of  the  positional  survey  data  was 
done  on  the  NPS  IBM  270/3033  AP  computer  system.  Graphic 
displays  were  constructed  using  the  Display  Integrated 
Software  System  and  Plotting  Language  (DISSPLA)  developed  by 
the  Integrated  Software  Systems  Corporation  (ISSCO) 

[Ref.  16].  All  computer  programs  involved  in  data 
processing  were  written  in  the  WATFIV  programming  language. 

Computations  were  made  in  an  X- Y  coordinate  system  based 
on  a  Modified  Transverse  Mercator  ( MTM)  projection.  A  MTK 
projection  is  essentially  the  same  as  a  Universal  Transverse 
Mercator  (UTM)  projection,  the  only  difference  being  that  in 
a  MTK  projection  a  central  meridian  is  picked  near  the 
survey  area  instead  of  being  fixed  at  a  particular  meridian 
[Ref.  17]. 

The  central  meridian,  controlling  latitude,  and  false 
easting  values  define  the  coordinate  system  used  for  compu¬ 
tations.  The  central  meridian  for  the  projection  was  chosen 
to  be  longitude  121°  52’  30"  W  which  is  approximately  the 
mean  longitude  of  the  survey  area.  The  controlling  lati¬ 
tude,  the  distance  ir  meters  from  the  eguator  to  a  reference 
latitude,  was  chosen  to  be  4,050,000  meters.  A  false 
easting  cf  5,000  meters  was  chosen  as  the  value  of  the 
X-coordinate  at  the  central  meridian. 

Three  shore  contrcl  stations  were  used  in  the  acquisi¬ 
tion  of  survey  data.  The  geodetic  positions  of  these 
stations  were  converted  to  the  X-Y  coordinate  system  (Table 
V)  using  program  UCOMPS,  which  is  a  hydrographic  utility 
package  available  to  students  at  NPS. 
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Range  information  was  recorded  using  a  Racal  Decca 
Trispccder  system,  a  microwave  system  commonly  used  for 
nearshore,  line-of- sight  survey  work.  On  October  28  arc 
November  30,  range-rarge  data  were  recorded  Ly  setting 
remote  units  over  stations  BEACH  LAE  and  MUSSEL.  Before  and 
after  the  survey,  the  ranging  system  was  calibrated  ever  the 
fixed  tase  line  USE  MCN  to  MUSSEL.  Daily  checks  in  the 
survey  area  were  made  to  determine  if  the  system  was  working 
properly.  This  was  accomplishea  by  maneuvering  the  survey 
vessel  tc  a  point  where  two  xnewn  navigational  ranges  inter¬ 
sected.  One  navigational  range  was  formed  by  stations 
MONTEEEY  AMERICAN  CAN  COMPANY  STACK  and  MONTEREY  RADIO 
STATION  K ME Y  MAST.  A  second  navigational  range  was  formed 
by  stations  MONTEREY  EAREOR  LIGHT  6  and  MONTEREY  BLUE 
LIGHTHOUSE. 

Track  control  for  range-azimuth  and  range-range  posi¬ 
tions  was  accomplished  by  steering  the  vessel  alcng  range 
arcs.  The  spacing  between  range  arcs  for  most  lines  was 
planned  to  be  40  meters.  Distance  between  positions  alcng  a 
sounding  line  averaged  approximately  200  meters.  The 
azimuth-azimuth  lines  were  controlled  by  steering  a  magnetic 
compass  heading. 

The  data  acquired  under  training  conditions  contained 
several  deficiencies  that  would  normally  not  be  tolerated. 
For  example,  the  quality  of  the  line  steering  was  generally 
poor;  the  vessel  wandered  off  the  arc  more  than  10  meters  in 
several  instances.  The  quality  of  the  sounding  lines  run 
using  azimuth-azimuth  control  was  extremely  deficient;  the 
position  plot  of  these  lines  show  a  jagged  path  by  the 
vessel.  Under  normal  hydrographic  procedures,  these  posi¬ 
tions  would  be  rejected.  Since  the  intent  of  this  study  is 
to  demonstrate  accuracy  analysis  techniques,  these  deficien¬ 
cies  prove  to  be  incc csequential ;  the  acquired  data  are 
adequate  to  demonstrate  the  concepts. 


Figure  3- 1  Hydrographic  Survey  Area 


ail  stations  are  of  third-order  or  better  and  are  published 
in  the  National  Geodetic  Survey  Data  Base. 

For  azimuth- azimuth  and  range-azimuth  positioning, 
azimuths  were  measured  with  a  Wild  T-2  theodolite.  Cn 
November  16,  range-azimuth  information  was  acquired  by 
locating  the  theodolite  over  station  MUSSEL  and  initialing 
cn  USE  MCN.  The  initial  direction  was  checked  by  sighting 
on  KM3Y  MAST.  Azimuth-azimuth  positions  were  acquired  cn 
November  23.  A  theodolite  was  set  over  USE  MON  and  an 
initial  direction  was  to  MUSSEL.  A  second  theodolite  was 
set  at  MCSSEL  using  USE  MON  for  the  initial  direction. 
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III.  EIPEEIHJNT  DESIGN  AND  IMPLEMENTATION 


The  coals  of  this  chapter  are  to  demonstrate  that  hydro- 
graphic  positioning  accuracy  can  be  classified  based  on  the 
radii  cf  SO  percent  confidence  circles  determined  by  using 
Hurt’s  method  and  to  show  that,  based  on  the  same  criteria, 
accuracy  predictions  can  be  made  for  survey  planning 
purposes. 

A.  DATA  ACQUISITION  IBOCEDUBES 

The  data  used  for  analysis  and  prediction  consisted  of 
range-range,  azimuth-azimuth  and  range-azimuth  survey  infor¬ 
mation.  The  data  were  acquired  by  Naval  Postgraduate  School 
(NPS )  students  in  a  Hydrographic  Sciences  course.  Although 
the  ccurse  was  structured  as  a  training  exercise,  the  data 
acquisition  procedures  utilized  were  nearly  identical  tc 
those  which  are  practiced  by  NOS. 

A  total  of  453  hydrographic  positions  were  recorded 
during  the  survey  of  a  nearshore  area  in  southern  Monterey 
Bay,  California.  Of  the  positions  used  for  analysis,  252 
were  range-range,  81  were  range-azimuth,  and  80  were 
azimuth-azimuth.  All  survey  information  was  recorded  tv 
hand  in  scunding  volumes.  The  vessel  used  was  a  36-foot 
Uniflite  with  a  fiberglass  hull  and  twin  engines.  The 
survey  was  conducted  cn  October  28,  November  16,  23,  and  30, 
1983.  Electronic  control  and  calibration  stations  used  for 
the  survey  included  CEE  MON  1978,  MUSSEL  1932,  BEACH  LAB 
1982,  MCNTEFEY  AMERICAN  CAN  CCMFANY  STACK  1932,  MCNTERE  Y 
EADIC  STATION  KMBY  MAST  1962,  MONTEREY  HARBOR  LIGHT  6  1978, 
and  WCNTEEEY  BLUE  LIGHTHOUSE  (Fig.  3.1).  With  the  exception 
of  MCNTEFEY  BLUE  LIGHTHOUSE,  which  is  a  low-order  position. 


Given  the  frequency  cf  1.6  MHz,  X  =  187.37  meters  frcm 

Equation  2.3.  The  lane  width  along  the  base  line  is  w'  =  w' 

g  r 

=  93.68  meters  from  Equation  2.5.  Using  the  law  of  cosines 

from  plate  geometry,  the  subtended  angles  a  and  a  are 

S  ^ 

32.47°  and  43.25°,  respectively .  The  angle  of  intersection 

cf  the  twc  hyperbolas  at  P  is  37.86°  from  Equation  2.7.  The 

lane  widths  at  P  are  w  =  254.  19  meters  and  w  =  335.  06 

r  g 

meters  from  Equation  2.6.  The  standard  errors  of  the  green 

(a  )  ana  red  {a  )  hyperbolas,  respectively  are  a  =  w  0  = 

i  :  l  g  base 

16.7  meters  and  o  =  w  a  =  12.7  meters.  These  standard 

2  r  base 

errors  are  in  a  linear  skewed  coordinate  system  and  must  be 

transformed  to  an  unccrrelated  rectangular  system.  Free 

Equations  2.18  and  2.19,  the  values  of  a  and  a  are  36.9 

a  b 

meters  and  12.7  meters,  respecti vely.  The  correlation  coef¬ 
ficient  in  the  correlated  rectangular  system  (P  b  )  is  then 
0.737  frcm  Equation  2.26.  The  semi-major  and  semi-minor 
axes  in  the  uncorrelated  rectangular  system  are  38.1  meters 
and  8.3  meters,  respectively,  from  Equations  2.27  and  2.26. 
The  eccentricity  is 

o 

c  =  ^  =  0.218 
x 

Table  IV  is  entered  with  the  values  of  P  =  0.9  and  c  = 

0.218.  The  value  for  K  is  found  to  be 

K  =  1.6602 

From  Equation  2.37,  the  radius  of  the  90  percent  probability 
circle  is  found  to  be 

H  =  63.3  meters 

The  probability  that  the  vessel's  position  will  be  within  a 
circle  of  63.3-meter  radius  centered  at  the  intersection  cf 
the  LCP's  is  90  percent. 
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o 


1.3  aeters 


a 

2 


c  =  -X  =  0.433 
a 

x 

Table  iv  is  er.tered  with  the  values  of  P  =  0.9  and  c  = 
0.433.  The  value  for  K  is  found  to  be 

7  =  1.7117 

Using  Equation  2.37,  the  radius  of  the  90  percent  prob¬ 
ability  circle  is  found  to  be 

P  =  5.14  meters 

The  prolal  ility  that  the  vessel's  position  will  be  within  a 
circle  of  5.14-aeter  radius  centered  at  the  intersection  of 
the  TCP's  is  90  percent. 


Example  3 

A  vessel  is  conducting  a  hydrographic  survey 
using  hyperbolic-hyperbolic  geometry.  The  hyperbolic  LC? 
generated  by  the  1.6-MHz  electronic  positioning  system  has  a 


star,  da  rl  error  of  0.95-lane  on  the  base  line. 


The  correla¬ 


tion  coefficient  (  p  )  betveer.  the  two  LCP's  is  known  tc  be 
0.4.  Compute  the  radius  of  the  90  percent  confidence  circle 
at  the  vessel’s  position. 

The  rectangular  plane  coordinates  of  the  master 
(T)  ,  two  slaves  (G  and  P)  ,  and  the  vessel's  position  (PI  are 


X  COORDINATE 


Y  COORDINATE 


172,679.  1 
209,  679.  1 
241, 738. 2 
223, 172.5 


t-2,540.  4 

99.540. 4 

21.325. 4 
169, 264. 2 


Exam pie  J 

A  vessel  is  conducting  a  hydrographic  survey 
using  range-range  geometry.  The  two  LOP's  generated  ty 

microwave  transmitters  have  standard  errors  of  a  =  3  meters 

1 

and  o  =4  meters.  The  angle  cf  intersection  8  at  the 

2 

vessel  is  30°.  Assuiie  the  LOP's  are  uncorrelated .  Compute 
the  probability  that  the  vessel’s  position  will  be  within  a 
circle  of  10-meter  radius  with  the  center  at  the  intersec¬ 
tion  cf  the  LOP's. 

Recalling  Equations  2.18  and  2.19,  the  values  of 

a  ana  a  are  found  tc  be  9.79  meters  and  6.14  meters, 
x  y 

respectively.  From  Equation  2.36 

o 

c  =  =  0.  633 

a 

x 

and  from  Equation  2.37,  with  R  =  10  meters, 

K  =  1.032 

Entering  Table  III  and  using  interpolated  values  for  c  and 
K,  the  probability  that  the  vessel's  position  will  be  within 
a  circle  of  10-meter  radius  centered  at  the  intersection  of 
the  ICP's  is 

P  =  53.2% 

Example  2 

A  vessel  is  conducting  a  hydrographic  survey 
using  range-azimuth  cecmetry.  The  range  LOP  generated  ly 
the  microwave  transmitter  has  a  standard  error  of  3  meters. 
The  azimuth  LOP  determined  by  theodolite  observation  has  a 
standard  error  of  1.3  meters  at  all  ranges.  Compute  the 
radius  of  the  90  percent  confidence  circle  at  the  vessel's 
posit io  r . 

In  the  range-azimuth  case  8  =  90°  and  the  LCP's 

are  uncorrelated.  Therefore, 

o  =  o  =  3.0  meters 
i  x 

and 
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TABLE  17 

Radii  cf  Circles  Given  c  and  P 


■“  \ 

0  0 

0.  1 

0.  2 

a  3 

0.  4 

0.  5 

0  ft 

0  7 

0.  8 

0  9 

1  0 

5000 

7  3  00 
■MHK) 

a  67449 
l  15035 

1  64485 

0  *3199 
l  15473 

1  04791 

ft.  70585 
l  16625 

1  65731 

0 

1 

74‘*'»3 

1*246 

67383 

0  '07.85 
l  231(H) 

1.  -j  99 18 

0  87042 

1  2'  ',34 

1  73708 

0  93365 

1  35.43 

1  7Q152 

(I  ‘)9<>2l 

1  42471 

1  S6253 

1  05769 
l  i023 1 

1.  94761 

I.  11807 

1.  5827 1 

2  04230 

1.  17741 

1.  >  651  ! 

2.  14597 

95.00 

j7:«) 

>900 

l  *5996 

2.  24 1  40 

Z  575».3 

1  9*1253 

2  24565 

2  57773 

1  97041 

2  25053 

2  5 a J?  7 

1 

98420 

'■*42  1 

2  00514 

2  2 8073 

2 

2  035  S*. 

2  .’.0707 

2  95257 

2  Oh i 30 
2  34581 

2  •■6533 

2  14*98 

2.  40156 

2  71515 

2  23029 

2  48494 

2.  76ub9 

2  331«0 

2.  58999 

2.  s9743 

2  44775 

2.  71620 

3.  J3455 

so 

*9 7  j 

2  '0703 
k  rd  j.»4 

I  29053 

2  808  w  3 
j 

j  2  Vi1  >6 

2  81432 

iJ10’,0 

1  *3 

3 

i 

8.1289 

0.'1888 

i  -4 -*9 

2  -  18.10 

3  » >214 

3  3!  "1  5 

2  85*94 

i  07144 

3  !3»*.4 

2  *8859 

3  .->671 

3  1*949 

2  93347 

1  15969 

3  39647 

J.  00431 

3  205*6 

3  4  56'** 

3  I  JC'73 

3.  31099 

3  55939 

3.  2.3525 
3  46164 
3  71692 

confidence  ellipse  is 

2 

A  =  Ka  o  tt 
e  x  y 

where  X  is  the  appropriate  probability  conversion 
(Table  I).  The  area  of  the  90  percent  confidence 

Ac  *  (2.39) 

where  E  is  given  by  Eguation  2.37.  For  a  condition  where 
o^  =  =  3  meters,  and  3  =  30°,  the  area  of  the  90  percent 

confidence  ellipse  is  261  square  meters,  while  the  area  of 
the  confidence  circle  is  587  sguare  meters.  For  both  stan¬ 
dard  errors  equaling  10  meters  and  8  =  30°,  the  90  percent 
confidence  ellipse  has  an  area  of  921  square  meters  and  the 
confidence  circle  has  an  area  of  2894  square  meters.  From 
an  operational  perspective,  the  difference  in  areas  between 
ellipses  and  circles  have  significant  implications  which 
will  be  discussed  in  Chapter  V. 

The  following  examples  are  presented  to  deacn- 
strate  methods  for  computing  the  parameters  cf  error 
ellipses  ard  confidence  circles  for  several  hydrographic 
positioning  geometries. 


(2.38) 

factor 
circle  is 


Ny  c 
K  N\ 

ao 

a  i 

0.  2 

a  3 

0.  4 

a  3  1 

1 

0  6 

07  j 

0.  8  1 

1 

0  9 

0.  1 

0790557 

0443987 

0242119 

.  0164176 

0123875 

.  009937"  ! 

0082940 

0071157  ! 

.  0062299  1 

.  0055400  ' 

0  2 

1585194 

.  13.39763 

0884533 

.  0628390 

0482413 

.  0390193 

0.127123 

0281415  1 

0246824  • 

0219757 

0  3 

2358228 

221.3804 

1739300 

.  1318281 

1039193 

.  085 1535  1 

07191C2 

0621386 

0540595  ■ 

0487639  ' 

0  4 

3108435 

3010228 

2635181 

.  2139034 

.  1742045 

1451S0S 

1237982 

107G237 

.  0950495  , 

08 50320 

0  5 

3829249 

.  5755584 

3481790 

.  3003001 

2532953 

2152S80 

1557448 

.  1C20829  J 

144394 l 

1290280  i 

0.  6 

4514938 

4457708 

4255605 

3840374 

3357384 

2914032  ! 

.  2548177 

. 2251114 

2009797  * 

18117*3  1 

0.  7 

5160727  , 

5115046 

4960683 

4633258 

4270862 

.  3699305  ■ 

3250302 

3935654 

2620373 

238:583 

0.  8 

.  5762892  ! 

5725957 

.  5004457 

.  53493S7 

4941382 

.  4474207  j 

4025028 

.  3627122  J 

3283453 

29*i07i*0 

0  9 

-  6318797  ! 

628.8721 

6191354 

.  5993140 

5051564 

52 3998  i 

.  4759375 

4333028  ! 

.  3953279  1 

3*' 20 135  i 

1.  0 

6826895  , 

.  6802325 

6723586 

.  6568242 

6291249 

5900953  i 

.  5461319 

.  5025790  ; 

4621421  | 

4257.553  ; 

1.  1 

7786079  1 

771.0597 

.  7707682 

.  7079681 

.  6859367 

6524480  j 

0116316 

.  5687497  ! 

.  5272402  j 

4887  373  : 

1.  2 

.  7698607  ; 

7682215 

7630305 

7532175 

.  7359558 

7079973  1 

0714269 

6306108  ! 

.  5893404 

.  54y*>736 

1  3 

SU63990  ! 

6050648 

8 00 8554 

.  7929968 

7793550 

.  7507365  ! 

7249073 

6873122  , 

6474394 

Go? ‘‘*22  i 

1  4 

8,3  S4 867  , 

8374049 

8340013 

8277048 

8169851 

.  7069288  1 

7720S89 

7353089 

7007000  . 

6023035 

1  5 

6663856 

. 8655127 

8627728 

8577362 

8493071 

.8350816  i 

81292S7 

.  7833962  j 

7459500 

7122546  i 

l  6 

8904014  1 

8897008 

8875060 

SS349U 

S70SG44 

8657559  ' 

8478393 

.  8220246 

.  7917194 

7574708  1 

1.  7 

9108691  , 

9103102 

9085619 

9053706 

9001746 

.  8915536  ! 

8773116 

8502471 

«291137 

7977882  , 

1  S 

9281394  ' 

9276964 

9203125 

9237989 

9197275 

9130080 

9019110 

8846024  , 

.  (5013238  ' 

.  8332175 

1  9 

9425669 

9422182 

9411299 

9391530 

9359855 

9308615  i 

<*222277 

9083609  j 

8*86731  ; 

8639149 

2  0 

.  9.544997  j 

9542272 

9533775 

9518415 

9*93815 

9454546  ] 

9388418 

.  9278709  | 

9:15702  , 

S 901495  ; 

2  1 

9642712 

.  9040598 

9634011 

9622127 

9603170 

.  9573305  j 

.  9522999 

9437668  i 

9305013  1 

9122714  1 

2.  2 

9721931  J 

9720304 

9715237 

.  9706109 

9691597 

.  9608845  ! 

9031017 

9565522 

9469380  . 

9306821  ! 

2  3 

9785518  : 

9784275 

9780408 

.  9773450 

9702419 

9745239  ' 

971*934 

96*)7306  i 

.  95*3739 

94 580® 5  ! 

2  4 

9836049  ’ 

. 9835108 

9832180 

.  982*918 

.  9813594 

.  9805703  ! 

9784661 

9747495  ■ 

9682698 

9580804 

2  5 

9875807  1 

. 9875100 

9872900 

.  9S68953 

9862720 

.  9853112 

9837569 

*810035  | 

9760522 

967  it  36 

.  9980*48 
.  9986257  ; 
9990332  I 
9993261 
.  9995347  ! 


9980542  ; 
99*6182 
9990279  , 
9993225  I 
9995323  I 


9897045  1 
9923483  I 
9943649  , 
.  9958878  | 
.  9970266 


.  9996505 
.  9907833  i  . 
999M12  I  . 
9998945  . 

.  9999305  i 


.  9099585 
.  9999732  ! 
.  999982S  1 
.  9999691 
.  9999932 


9994218  j  . 
. 9990102  . 

9997396 
.  9f»‘i827G 
9998570  j  . 

.  9999266  ,  . 
.  9999527  '  . 
.  9999698  ;  . 
. 9999809  j  . 
.  9999681  |  . 


9999955  |  . 
9999973  . 

99099*3  | 
9999990  ;  . 
9999994  |  . 
I 

9999996 
9999908  ! 

.  9999999  I  . 


.  9999951  . 

. 9999970 
.  9999982  1 
.  9099989  1 
.  0900993  ! 


9999996  i  . 
'.'999098  !  . 
J900999 


5.  4 

9909999 

0999999 

9909999  |  9999909  i  9909999  !  9909999  1 

. 9999990 

9099999  i 

9999999  , 

999099 s 

9' '9999  5 

j.  5 

1.  OOOOOOO 

1.  ooooooo 

1  OOOOOUO  .1.  UU00000  !  1.  UOOOOOO  l.  0000000 

J  0000000 

91I99999  1 

.  9999999  j 

9999999 

9999997 

5.  6 

1  1 

1.  <KK)0000  1 

OOOOOOO  ■ 

.  9999999  1 

9J9999S 

5  7 

I 

1 

! 

1  UOOOOW) 

999.1909 

5  8 

l 

1 

1  iXXKhiOO 

5  9 

'  1 

i 

j 

6  0 

!  1 

J 

i 

i 

i 

_ !_ 

I 


b.  Circles  cf  Equivalent  Probability 


Burt  [Eef.  3]  presents  a  method  for  translating 
ellipses  cf  equivalent  probability  into  circles  of  equiva¬ 
lent  probability.  Tc  utilize  this  method,  it  is  first 
necessary  to  compute  the  eccentricity  cf  the  error  ellipse, 
c,  by  the  equation 


c 


a 


JL 

a 

x 


(2.36) 


where  a  >  a  . 

x  y 

Harter  [Bef.  15]  compiled  Tables  III  and  IV 
which  are  taken  from  Bowditch  £  Ref .  10,  p.  1215]-  Harter's 

data  are  given  in  terirs  of  the  eccen tricity,  c,  a  parameter, 
X,  and  a  probability,  P.  The  parameter,  X,  when  multiplied 
by  ax  gives  the  value  of  the  radius,  R,  of  the  circle  of  the 
corresponding  probability  shown  in  Table  III.  That  is, 

"-*°x  12.37) 

The  probability  of  a  point  falling  inside  a 
circle  cf  specified  radius  can  be  computed  by  entering  Table 
III  with  c  and  K  as  arguments.  Given  a  fixed  probability,  K 
is  determined  by  entering  Table  IV  using  c  and  P  as  argu¬ 
ments.  The  radius  of  the  probability  circle  is  then 
computed  using  Eguaticn  2.37. 

Using  confidence  ellipses  has  certain  advantages 
ever  confidence  circles  of  equal  probability.  First,  the 
directional  nature  of  the  true  error  distribution  is  net 
represented  in  the  confidence  circle  method  even  though  both 
methods  give  an  accurate  measure  of  confidence.  Second,  the 
area  of  the  confidence  ellipse  is  always  less  than  or  equal 
to  the  area  of  the  confidence  circle.  The  area  of  a 
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i.  C.5  mm  at  the  scale  of  the  survey  for  scales  of 
1:20,000  and  smaller, 

ii.  1.0  mm  at  the  scale  of  the  survey  for  1:10,000 
scaie  surveys,  or 

iii.  1.5  mm  at  the  scale  of  the  survey  for  scales  of 
1:5,000  and  larger. 

The  major  advantage  of  using  d  as  a  precision 

rms 

index  is  its  ease  of  computation.  Some  hydrographers  draw 

analogy  between  the  varying  probability  associated  with  ore 

d  (€3.2  percent  to  68.3  percent)  and  the  fixed  prct- 
rms 

ability  associated  with  a  one-dimensional  standard  error 

(68.3  percent).  In  fact,  d  has  very  little  statistical 

rms 

meaning.  The  obvious  problem  with  using  dras  as  a  precision 
index  is  the  varying  probability  associated  with  the  error 
circle.  Fcr  this  reason  Greenwalt  [Eef.  12,  p.  31]  recom¬ 
mends  against  its  use. 


TABLE  II 
Probabilities  Associated  iith  d 


rms 


PROBABILITY 


0 

y 

a 

X 

LENGTH  OF 

1  ^s 

1  *rms 

0.0 

1.0 

1.  000 

0.683 

0.  1 

1.0 

1.  005 

0.632 

0.2 

1.0 

1.020 

0.682 

0.3 

1 .  0 

1.  042 

0.6  76 

0.4 

1.0 

1.  077 

0.671 

0.5 

1.0 

1.  118 

0.662 

0.6 

1.0 

1.  166 

0.650 

0.7 

1.0 

1.  220 

0.64  1 

0.8 

1.0 

1.  280 

0.635 

0.9 

1.0 

1.  345 

0.632 

1.0 

1.0 

1.  4  14 

0.632 

2  d 


rms 


0.  954 
0.955 
0.957 
0.961 
0.  966 
0.  969 
0.973 
0.  977 
0.  980 
0.  981 
0.  982 


An  error  circle  with  a  radius  of  one  d  car.  be 

rms 

constructed  about  the  intersecting  LOP's  (Fig.  2.9).  Two 

d  is  the  radius  cf  the  error  circle  obtained  using  two 
rtns 

times  the  values  of  a  and  a  in  Equation  2.31.  For  an 

x  y 

elliptical  error  distribution,  the  probability  associated 

with  a  specific  value  of  d  varies  as  a  function  of  the 

rms 

eccentricity  of  the  error  ellipse  (Table  II).  The  prob¬ 
ability  associated  with  one  dras  varies  from  63.2  percent  to 
68.3  percent,  while  the  probability  associated  with  twc  d 

rms 

varies  between  95.4  percent  and  98.2  percent. 


distribution  can  be  substituted  for  the  elliptical  distribu¬ 
tion.  This  substitution  can  be  satisfactory  for  error 
analysis  within  certain  °y'/°x  ratios.  However,  when  this 
ratio  is  small  the  distortion  introduced  by  the  circular 
distribution  may  become  misleading. 

a.  Hoot  Mean  Square  Error 

The  terms  radial  error,  root  mean  square  error, 

and  d  are  identical  in  meaning  when  applied  to  two- 

rms  3  rr 

dimensional  errors  [  Hef .  10,  p.  1229  ].  The  term  d  is 
defined  as  the  square  root  of  the  sun  of  the  squares  of  the 
standard  errors  along  the  major  and  minor  axes  of  the  error 
ellipse.  That  is 


d  «  \/cr  +  a 

rms  x  y 

(2. 

where  a ^  and  are  given  by  Equations  2.18  and  2.19. 

A  more  direct  form  of  2.31  is  given  by  [Hef.  2,  p.  54] 


(2.21) 


sins  V  i 


a2  +  a2 


(2.32) 


for  uncorrelated  LOP’s.  For  range-range  and  azimuth-azimuth 
positioning,  with  CTj  =  °2  =  0  •  Equation  2.32  reduces  tc 


d-  -  &  a 

rms  sins 


(2.33) 


For  range-azimuth  positioning,  B  =  90°  and  Equation  2.32 
becomes 


(a1  +  a2 

l  2 


(2.34) 


The  mere  general  form  of  Equation  2.32  for  both  correlated 
and  uncorrelated  LOP’s  [Hef.  2,  p.  59]  is 


drms  *  ITHs  V°!+  a 2  *  2p12°,a2COs8 


(2.25) 


where  0  is  the  correlation  coefficient. 


r 


V 

k 


a 

a 


used  to  compute  I-Y  coordinates  of  all  positions.  Based  cn 
geometric  relationships  discussed  earlier,  angles  of  inter¬ 
section  of  the  ICP's  were  then  computed  for  range-range  ar.d 
azimuth-azimuth  points.  The  angles  of  intersection  for  all 
range-azimuth  positions  are  90°. 

The  range-range  ard  azimuth-azimuth  data  were  then 
passed  to  WATFIV  subroutine  PRCB  (Appendix  A).  As  input 
parameters,  the  subroutine  accepts  two  standard  errors  cf 
the  ICP’s  and  the  corresponding  angle  of  intersection.  The 
output  parameters  include  the  semi-major  and  semi-minor  axes 
cf  the  90  percent  corfidence  ellipse,  the  radius  of  the  90 
percent  confidence  circle,  and  the  areas  covered  by  both 
figures. 

Subroutine  PEOB  uses  a  linear  approximation  to  determine 
the  value  of  the  function  K  for  varying  values  of  the  eccen¬ 
tricity,  c,  in  Burt's  method.  A  linear  interpolation  was 
performed  by  first  taking  the  eleven  discrete  values  of  c 
and  K  for  a  probability  of  90  percent  from  Table  IV  and  then 
constructing  a  series  of  relationships  for  K  as  a  function 
cf  c  (Table  VI) . 

Values  of  the  radii  of  90  percent  confidence  circles  for 
range-range  data  were  plotted  at  their  respective  positions 
(Fig.  4.1).  The  arcs  of  circles  connecting  the  two  control 
stations  BEACH  LAB  and  MUSSEL  represent  lines  of  constant 
intersection  angle  (30°).  Of  the  range-range  data  set, 
position  948  (Appendix  3) — coordinates  X  =  4119.01,  Y  = 

4735. C7 — was  found  to  have  the  smallest  radius  (strongest 
position)  of  6.4  meters  and  an  angle  of  intersection  cf 
90.2°.  Fosition  1 37--coordi na tes  X  =  3345.86,  Y  = 

3873. 34--represents  the  weakest  position  with  radius  value 
of  15.3  aeters  and  an  angle  of  intersection  of  26.7°. 

The  positional  accuracy  degrades  rapidly  as  the  inter¬ 
section  angle  approaches  30°;  the  30°  arc  represents  a  line 
of  constant  13.7  meter  radius.  Within  400  meters  of  the  50° 
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1 

TABLE  VI 

Linear  Approximations  for  K  as  a  Function  of  c 

Interval  of  c  Linear  Interpolation 

~  ~  TUncf ion  for  K 


0.0 

- 

0.  1  | 

1  K 

= 

.0306c 

♦ 

1  .  64485 

0.1 

- 

0.2 

K 

.0940c 

+ 

1.63851 

0.2 

- 

0.  3 

K 

= 

.  1652c 

•fr 

1.62427 

0.3 

— 

0.4 

K 

= 

.  2535c 

•f 

1 . 59773 

0.4 

- 

0.5 

K 

= 

.  3790c 

♦ 

1 . 54753 

0.5 

- 

0.6 

K 

= 

.  5444c 

♦ 

1 .46433 

0.6 

- 

0.7 

K 

= 

.  710  1c 

♦ 

1.36546 

0.7 

- 

0.8 

K 

.  8508c 

♦ 

1.26697 

0.3 

- 

0.9 

K 

.  94  7  5c 

♦ 

1  .  1896  1 

C.  9 

- 

1.0 

K 

= 

1. 036  1c 

1 .10937 

intersection  arc,  the  radius  varies  between  8  and  15  meters.  • 

The  radii  values  charge  slowly  in  the  vicinity  of  the 
minimum  value  of  6.4  meters  which  corresponds  to  an  angle  of 
intersection  of  90°. 

The  radii  of  90  percent  confidence  circles  associated  • 

with  the  azimuth-azimuth  positions  acquired  using  control 
stations  USE  MON  and  PUSS EL  were  also  plotted  at  their 
respective  positions  (Fig.  4.2).  The  standard  errors  cf  the 
LOP’s  are  assumed  to  te  1.3  meters;  the  resulting  improved  • 

accuracy  is  evident.  The  maximum  value  of  the  90  percent 
confidence  circle  radii  is  8.7  meters  at  position 
637--coordinates  X  =  4327.25,  Y  =  28 18. 39--w hich  corresponds 
to  an  angle  of  intersection  of  159.8°  (or  in  terms  of  the  9 

supplement,  20.2°).  Eosition  682 — coordinates  X  =  4611.20, 

Y  =  4421.29 — represents  the  strongest  position  recorded 

during  the  survey  with  a  90  percent  confidence  circle  radius 

of  2.8  meters  and  an  angle  of  intersection  of  91.0°.  ^ 

Again,  the  rapid  degradation  of  accuracy  is  noted 
approaching  S  =  150°.  The  arc  of  the  150°  intersection 
angle  represents  a  constant  radius  of  5.9  meters.  Discrete 
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Figure  4. 1  Range-Range  Accuracy  Analysis 
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values  along  the  arc  confirm  this  qualitatively.  A  large 
area  of  strong  positional  accuracy  surrounds  the  area  where 
8  =  SO0.  Numerous  values  of  2.8  meters  are  present  near 

the  top  cf  the  plot. 

Using  the  assumptions  stated  at  the  beginning  of  this 
section,  the  values  for  all  radii  of  90  percent  confidence 
circles  for  range-azimuth  positions  are  5.  1  meters.  This 
computation  was  carried  out  in  Example  2  of  Chapter  II. 

Since  this  case  is  trivial,  the  data  are  not  displayed 
graphically. 

Positioning  data  were  also  classified  cased  on  the 
parameters  of  the  90  percent  confidence  ellipse.  WATFIV 
program  E1LIP  (Appendix  C)  was  used  to  generate  the  parame¬ 
ters  cf  the  90  percent  confidence  ellipse  for  range-range, 
azimuth-azimuth,  and  range-azi uuth  positioning  data.  The 
program  was  initialized  by  entering  the  coordinates  cf  the 
control  stations  and  standard  errors  of  the  LOP’s.  The  fix 
number,  hydrographic  position  coordinates,  and  angle  of 
intersection  were  then  read  in  from  a  data  file.  Subroutine 

PEOB  was  called  to  compute  values  for  Ka  and  K o  . 

x  y 

The  angle  of  orientation  of  the  major  axis  of  the 
ellipse,  measured  clockwise  from  north,  was  then  computed, 
for  range-range  and  azimuth-azimuth  positions,  the  LCP 
generated  from  the  left  control  station  was  used  as  the  base 
LOP.  For  range-azimuth  positions,  the  LOP  formed  by  the 
theodolite  was  used  as  the  base  LOP.  First,  the  orientation 
of  the  base  LOP  in  the  coordinate  system  was  determined. 

The  orientation  of  the  major  axis  of  the  error  ellipse  rela¬ 
tive  to  the  base  LOP  (9)  was  then  computed  using  Eguaticn 
2.29.  By  adding  or  subtracting  9  to  the  orientation  of  the 
base  LOP,  the  orientation  of  the  major  axis  of  the  error 
ellipse  in  the  coordinate  system  was  determined.  This  angle 
takes  on  values  from  0°  to  180°.  Appendix  D  consists  cf  the 
confidence  ellipse  classification  scheme  for  range-range. 
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azimut h-azimuth  and  range-azimuth  data.  Forty  positions  for 
each  positioning  geometry  are  listed  for  comparison  to  the 
classification  scheme  presented  in  Appendix  3. 

Appendix  B  lists  the  data  fcy  position  number,  X-Y  coor¬ 
dinate,  angle  of  intersection,  and  radius  of  the  90  percent 
confidence  circle.  Appendix  D  lists  the  data  by  position 

number,  X-Y  coordinate,  angle  of  intersection,  Ka  ,  Ka  ,  and 

x  y 

angle  of  orientation  for  the  90  percent  confidence  ellipse. 
These  appendices  are  similar  to  hydrographic  survey  data 
bases  and  demonstrate  accuracy  classification  schemes  based 
on  the  two  criteria. 

C.  ACCURACY  PREEICT ICNS 

The  overall  positional  accuracy  of  a  survey  can  be 
controlled  by  computing  accuracy  values  before  data  acquisi¬ 
tion  is  begun.  For  example,  if  the  hydrographer  is  using 
radii  of  90  percent  confidence  circles  as  an  accuracy 
criterion,  the  minimum  allowable  angle  of  intersection  for 
two  ICF's  can  be  computed  for  meeting  specif ications.  The 
nature  of  the  survey  area  may  allow  the  flexibility  to 
change  system  geometry  to  maximize  accuracy  at  a  specific 
location  or  to  maximize  the  area  covered  with  a  giver,  accu¬ 
racy.  By  making  accuracy  computations  before  acguiring 
data,  the  hydrographer  may  also  have  the  option  of  deciding 
what  type  of  positioning  system  is  to  be  used  to  meet  accu¬ 
racy  requirements. 

The  construction  cf  reliability  contours  is  one  method 

to  display  the  expected  positional  accuracy.  Reliability 

contours,  lines  cf  constant  repeatable  accuracy  which  are 

functions  of  the  system  geometry  and  standard  errors  cf  the 

positioning  equipment,  can  be  constructed  about  shore 

stations  using  the  radii  of  90  percent  confidence  circles 

criterion  or  the  less  desirable  d„.  value. 

rras 


59 


Consider  the  equations  that  have  been  developed  in 
Chapter  II  for  the  determination  of  radii  of  90  percent 
confidence  circles  using  Burt's  method.  For  uncorrelated 
LOP's  in  a  range-ranee  or  azimuth-azimuth  system,  the 
repeatable  accuracy  cf  a  hydrographic  position  is  a  function 
only  cf  the  angle  of  intersection,  assuming  the  standard 
errors  of  the  LOP's  are  constant  throughout  the  survey  area. 
The  locus  of  points  which  define  a  constant  angle  of  inter¬ 
section  for  two  LOP's  in  a  range-range  or  azimuth-azimuth 
system  is  a  circle  which  passes  through  both  control 
stations.  Given  the  coordinates  of  the  two  control 
stations,  the  equations  of  these  circles  can  be  determined. 

Construction  of  reliability  contours  involves  several 
simple  trigonometric  relationships  (Fig.  4.3).  Let  LR  be 
the  line  connecting  the  two  shore  control  stations  L  and  E 
in  a  range-range  system.  The  length  of  line  LR  is  b.  The 
circle  through  both  stations  defines  a  line  of  constant 
intersection  angle  for  two  LOP's.  The  radius  of  the  circle 
is  r.  The  distance  e  is  measured  along  the  perpendicular 
bisector  of  the  line  IE  to  the  center  of  the  circle  at 
point  0(h,k)  and  is  given  by 

b 

e  *  2tan8  (4.1) 

Knowing  e  and  the  radius  r,  the  coordinates  of  point  0  can 
be  computed.  The  equation  of  the  circle  is  then 

r1  *  (x  -  h)2  +  (y  -  k)2  (4>2) 

These  two  equations  were  used  to  generate  reliability 
contours  for  display  on  a  computer  graphics  terminal.  Using 
Eurt's  method,  the  angles  of  intersection  of  two  LOP's  were 
computed  for  discrete  values  of  radii  of  90  percent  confi¬ 
dence  circles.  Reliability  contours  about  stations  EEACH 
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Figure  4.3  Construction  of  a  Seliability  curve 


IA3  and  MUSSEL  for  a  range-range  system  (a  =a  =  3  meters) 

1  2 

were  constructed  (Fig.  4.4).  Using  Equation  ‘4.2,  X-Y  coor¬ 
dinates  were  generated  for  points  laying  on  different  reli¬ 
ability  circles.  A  curve-fitting  subroutine  in  the  EISEPIA 
library  was  used  to  generate  the  circles  through  the 
computed  points.  The  13-meter  accuracy  contour  corresponds 
to  an  angle  of  intersection  of  31.6°,  while  the  7-meter 
accuracy  contour  corresponds  to  an  angle  of  intersection  of 
67.9°.  The  best  achievable  accuracy  of  the  system  is  6.4 
meters  at  90°. 

Fcr  comparison  purposes,  reliability  contours  were 

constructed  about  BEACH  LAB  and  MUSSEL  for  azimuth-azimuth 

geometry  (  a  =a  -  1.3  meters).  The  increased  accuracy  of 
l  2 

this  configuration  is  evident  (Fig.  4.5).  The  3-meter 


contour  corresponds  tc  an  angle  of  intersection  of  69. 4C 
while  the  6-meter  contour  corresponds  to  an  angle  cf  inter¬ 
section  of  29.6°.  The  best  achievable  accuracy  at  an  inter¬ 
section  angle  of  90°  is  2.8  meters. 

A  second  scheme  was  used  tc  display  accuracy  predictions 
for  the  two  positioning  methods.  Given  the  coordinates  cf 
E2ACH  IAE  and  MUSSEL,  a  series  of  discrete  points  spaced  800 
meters  apart,  were  generated  throughout  the  survey  area. 

The  values  for  the  radii  of  90  percent  confidence  circles 
were  theD  computed  at  each  point  with  the  use  of  subroutine 
ER03 .  Figures  4.6  and  4.7  illustrate  this  prediction 
scheme.  These  figures  present  the  same  information  as 
Figures  4.4  and  4.5  in  a  different  manner.  The  30°  angle  of 
intersection  contour  is  shown  on  both  figures. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


A.  ACCURACY  SPECIFICATIONS 

Interpretation  or  the  1932  IKC  positioning  standards  in 
terms  cf  90  percent  confidence  circles  yields  some  inter¬ 
esting  results  with  respect  to  present  day  survey  practices. 
For  example,  for  a  1  :  10,000-scale  hydrographic  survey,  NOS 
usually  uses  microwave  positioning  systems  in  a  range-range 
mode,  and  assumes  a  standard  error  of  3  meters  for  each  LCP. 
Surveys  are  frequently  conducted  between  the  30°  tc  150° 
angle  of  intersection  limits.  Using  the  90  percent  confi¬ 
dence  circle  critericr,  the  radius  of  the  circle  should  net 
exceed  10  meters.  However,  the  radius  value  for  8  =  30°  and 

150°  is  13.7  meters.  The  values  of  K  a  and  Kcr  for  the  90 

x  y 

percent  confidence  ellipse  are  17.6  and  4.7  meters,  respec¬ 
tively.  To  meet  the  SO  percent  criterion  for  a  1:10, 3  30- 
scale  survey,  the  8  limits  should  be  42°  to  133°. 

Azimuth-azimuth  positioning  is  accurate  enough  for 
1; 5 , 0 0 C-scale  surveys,  using  8  limits  of  35°  to  145°, 
assuming  a  standard  error  of  1.3  meters  for  each  LCP.  with 
the  standard  error  assumptions  used  for  range-azimuth,  the 
90  percent  radius  is  5.1  meters  for  all  positions.  Given 
the  uncertainties  of  the  standard  error  figures,  it  is 
rational  tc  assume  that  range-azimuth  positions  can  meet  tne 
5-meter  accuracy  standard  for  1 : 5 , 000-scale  surveys.  In 
fact,  range-azimuth  positional  accuracy  can  exceed  azimuth- 
azimuth  accuracy  when  the  later's  S  is  less  tnan  35°.  For  a 
3-meter  a  range-range  configuration,  it  is  impossible  tc 
meet  1:5,000  specifications  with  any  8. 

As  a  general  guideline,  the  30°  to  150°  angle  of  inter¬ 
section  limit  is  a  geed  rule  tc  use  for  uncorrelated  LCE's. 
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CLASSIFIED  AZIMUTH-AZIMUTH  POSITIONS 


Control 

Stations: 

USE  MCN  1978 

and  MUSSEL  1932 

Stan  card 

Error  Used 

ir  Compu ta tions:  1.3  meters 

Fix 

X 

Y 

Anqle  of  Radius 

No .  Coordinate 

Coordinate 

Intersection  90^  ci 

615 

4449.26 

27  1 1.40 

156.  9 

7.6 

620 

4437. 86 

2807. 29 

153.0 

6.5 

621 

4427. 4  1 

2897.  17 

149.6 

5.  9 

622 

4419.65 

2989.36 

146.2 

5.  3 

623 

4410.  42 

3072.  20 

143.4 

4.  9 

624 

4390. 60 

3154.04 

141.3 

4.7 

625 

4376. 52 

3234.  24 

138.9 

4.  4 

626 

4364.  3  1 

33 19. 82 

136.2 

4.2 

627 

4348.  7  1 

3411.24 

133.  5 

4.  0 

628 

4334.31 

3502.  51 

130.7 

3.3 

631 

4339.97 

3331. 55 

135.  1 

4.  1 

632 

4337. 95 

3290. 68 

138.7 

4.4 

633 

4327.70 

3199.92 

142.  8 

4.3 

634 

4322. 53 

3107.83 

146.  9 

5.4 

635 

4323.  1  8 

3012.49 

151.0 

6.1 

636 

4324.42 

2916.00 

155.  3 

7.  1 

637 

4327.  25 

2818.39 

159.  8 

8.7 

638 

4394.77 

2806. 86 

156.0 

7.3 

639 

4336. 58 

2903.  09 

152.  0 

6.3 

640 

4377. 29 

2993.63 

148.4 

5.6 

64  1 

4367. 00 

3090. 26 

145.  1 

5.  1 

642 

4355.  55 

3187. 06 

141.9 

4.7 

643 

4345.  97 

3285.  88 

133.  5 

4.4 

644 

4256.  90 

3516. 02 

133.6 

4. 0 

645 

4260.68 

3416. 65 

137.  4 

4.  3 

646 

4264. 77 

3321. 01 

14  1.1 

4.7 

647 

4283.59 

3208. 82 

144.  8 

5.  1 

648 

4293. 24 

3127.  36 

147.  7 

5.5 

649 

4300.30 

3024.  70 

151.8 

6.3 

650 

4345.  53 

3  145.  58 

144.  1 

5.  1 

651 

4370.38 

3236.  09 

139.  1 

4.5 

652 

4398.77 

3327. 77 

134.2 

4.0 

653 

4411.48 

3421. 46 

130. 2 

3.3 

654 

4438.  30 

3506. 23 

125.9 

3.5 

655 

4470. 97 

3591.21 

121.0 

3.3 

656 

4502. 73 

3677. 50 

117.4 

3.2 

657 

4514.38 

3767.  26 

114. 0 

3.  1 

658 

4512. 00 

3860. 1 3 

111.1 

3.0 

65  S 

4520.  10 

3948.53 

107.9 

2.9 

660 

4494.79 

3049.  75 

139.2 

4.5 

66  1 

4487. 96 

3  14  4.  60 

136.  2 

4.2 

6  6  2 

4477. 62 

3243. 86 

133.  2 

4.0 

663 

4462. 90 

3372. 74 

129.4 

3.  7 

664 

4453. 74 

3469. 99 

126.5 

3.5 

665 

4440. 56 

3564. 89 

123.  8 

3.4 

666 

4465-  02 

3652. 95 

119.7 

3.  3 

667 

4507.24 

3743.  29 

115.0 

3.  1 

668 

4577.55 

3595.  27 

116.9 

3.2 

669 

4569.  50 

363  1.0  2 

114.6 

3.  1 

670 

4563. 68 

3776.77 

111.9 

3.0 

67  1 

4564.53 

3872.  77 

108.  3 

3.0 

672 

4563.90 

3965.  68 

106.0 

2.9 

673 

4560.  4  1 

4057.  10 

103. 2 

2.  9 

674 

4553.  49 

4 147. 06 

100. 6 

2.8 

675 

4556.  7  1 

4239. 64 

97.  7 

2.  3 

677 

4571.60 

4416.22 

92.0 

2-8 

678 

4576.20 

4504.  92 

89.3 

2.8 

679 

4582. 23 

4597.03 

36.5 

2.3 

690 

4604.35 

4631.  33 

85.  1 

2.  8 

68  1 

4613. 39 

4527.  1  7 

87.9 

2.  8 

RANGE-RANGE  ACCURACIES  (CONTINUED) 


Fix 

X 

Y 

Angie  of 

Radius  of 

No. 

Coordinate 

Coordinate 

Intersec  tion 

90%  circle 

806 

3630. 54 

4843. 53 

92.3 

6.4 

307 

3877. 77 

4734. 1 8 

88.3 

6.4 

808 

3907.61 

46  18.  1  8 

33.4 

6.  5 

809 

3932.77 

4495.  95 

73.  1 

6.6 

8  10 

3948.  19 

4376. 06 

72.  8 

6.8 

811 

3956.79 

4251. 77 

67.0 

7.  0 

8  12 

3953.  8  1 

4125.89 

60.9 

7.  5 

813 

3938.32 

3995.  79 

54.  1 

8.  1 

814 

3909. 94 

3869. 78 

46.  9 

9.  1 

815 

3864.71 

3744. 84 

33.9 

10.  7 

816 

3816. 98 

3619-07 

30.6 

13.  4 

817 

3711.62 

5131.00 

104.  5 

6.7 

318 

3788. 16 

5014.61 

99.7 

6.5 

819 

3851. 30 

4889.79 

94.  7 

6.  5 

820 

3901.  13 

4761. 08 

39.6 

6.4 

821 

3961.73 

4590. 52 

82.7 

6.5 

822 

3977.  10 

4485.23 

78.2 

6.6 

323 

3993.32 

4  1  18.  3  1 

61.5 

7.4 

824 

3996. 85 

4193.  81 

65.2 

7.2 

825 

3981.72 

4048. 49 

57.  9 

7.7 

826 

3956.  88 

3905.60 

50.2 

8.  6 

827 

3916. 26 

3758.91 

41.  5 

10.1 

828 

3826. 04 

3553. 55 

27.6 

14.  a 

829 

3961.88 

4691.  6  1 

87.  1 

6.4 

830 

3807. 91 

4827.37 

92.  4 

6.4 

831 

3841. 30 

4953.  23 

97.3 

6.5 

832 

3789. 42 

5085.46 

102.4 

6.6 

333 

3772.12 

5177.05 

105.  a 

6.7 

834 

3843.  76 

5069- 86 

101.7 

6.6 

835 

3906.65 

4961.  46 

97.6 

6.  5 

836 

3959.36 

4844. 18 

93.3 

6.4 

837 

4002.25 

4720.39 

88.6 

6.4 

838 

4057. 98 

4659.  81 

86.7 

6.  4 

339 

3992.26 

4828.76 

92.8 

6.4 

84  C 

3954. 08 

4915.  16 

96.0 

6.  5 

841 

3905. 75 

5046.  63 

100.  8 

6.6 

842 

3822. 81 

5167. 39 

105.  2 

6.7 

843 

3844.43 

5129.  36 

103.  8 

6.6 

844 

3858-  55 

5033.  27 

100.3 

6.6 

845 

4016. 73 

4919.44 

96.  4 

6.  5 

846 

4061. 27 

4797.1 3 

92.  1 

6.  4 

847 

4159. 72 

4597.99 

35.5 

6.8 

848 

4 1 19. 0  1 

4735.07 

90.2 

6.4 

849 

4077. 07 

4872.  84 

95.  0 

6.5 

850 

4021. 39 

4990. 89 

99.0 

6.5 

85  1 

3961. 68 

5118.  03 

103.  3 

6.6 

852 

3992. 88 

5139.94 

104.0 

6.7 

853 

4056.5  1 

5026.  1  5 

100.  3 

6.6 

854 

4121. 60 

4874.08 

95.4 

6.  5 

855 

4154. 57 

4779.44 

92.2 

6.  4 

856 

4196. 42 

4617. 75 

86.7 

6.4 

RANGE-RANGE  ACCURACIES  (CONTINUED) 


Fix 

X 

Y 

Angle  of 

Radius  cf 

No. 

Coordinate 

Coordinate 

Intersect  ion 

90%  Circle 

742 

3570.  50 

4836. 50 

92.8 

6.4 

743 

3618. 79 

4723.  26 

87.0 

o.4 

744 

3662. 90 

4602. 44 

80.9 

6.  5 

745 

3691. 10 

4478.  20 

74.4 

6.7 

746 

3715.24 

4341.39 

67.  1 

7.0 

747 

3720. 11 

4205. 94 

59.5 

7.5 

748 

3703. 21 

4071. 17 

51.4 

8.5 

749 

3678.03 

3939.55 

43.  0 

9.8 

7  5  C 

3636. 04 

3801.  29 

33.7 

1  2.  3 

751 

3630. 00 

3715. 02 

28.  3 

14.5 

752 

3678. 77 

3812.  88 

35.7 

11.6 

753 

3715.73 

3925. 56 

43.  4 

9.7 

754 

3731.4  8 

4050. 62 

51.0 

8.5 

755 

3756.38 

4168.  71 

58.3 

7.7 

756 

3754. 68 

4300. 89 

65.5 

7.  1 

757 

3743. 02 

4427. 92 

72.2 

6.8 

758 

3710. 29 

4554. 17 

78.6 

6.6 

759 

3680.62 

4673. 88 

84. 6 

6.5 

760 

3627.42 

4792. 5 1 

90.5 

6.4 

76  1 

3572.80 

4899. 46 

95.  8 

6.  5 

762 

3627. 15 

4891.25 

95.1 

6.5 

763 

3683.33 

4778.34 

89.  7 

6.4 

764 

3725.77 

4656.0  1 

83.9 

6.  5 

765 

3757. 98 

4525.66 

77.5 

6.  6 

766 

3781.56 

4403.  1  3 

71.5 

6.8 

767 

3798.56 

4272.  52 

64.8 

7.2 

768 

3787. 95 

4134.  55 

57.  1 

7.3 

769 

3772.73 

4001.69 

49.4 

8.7 

770 

3738. 73 

3870. 33 

41.  1 

10.2 

771 

3695.62 

3745.  45 

32.  6 

12.7 

772 

37 12- 52 

3708.43 

31.2 

13.2 

773 

3777.  45 

3829.  70 

40.  2 

70.4 

774 

3803. 43 

3957. 03 

48.0 

9.0 

775 

3818.28 

4087. 24 

55.  4 

8.0 

776 

3831. 33 

4211.  74 

62.3 

7.4 

777 

3823.70 

4338.97 

68 ;  8 

6.9 

778 

3813.4  1 

4476. 12 

75.6 

6.7 

779 

3777. 13 

4614.24 

82.  1 

6.5 

780 

3739.58 

4736.39 

87.8 

6.4 

781 

3691.63 

4868. 16 

93.9 

6.  5 

782 

3706. 65 

4918.  98 

96.0 

6.  5 

783 

3749.46 

4806. 82 

91.  1 

6.4 

784 

3799. 89 

4705.  4  1 

86.6 

6.4 

785 

3829.30 

4581. 1 4 

81.0 

6.5 

786 

3857. 30 

4455. 00 

75.  1 

6.7 

787 

3868.49 

4327.  66 

69.0 

7.0 

788 

3872. 94 

4198. 56 

62.5 

7.3 

789 

3862.59 

4062. 95 

55.  3 

8.0 

790 

3841. 17 

3928.97 

47.  7 

9.0 

791 

3804.62 

3797.  14 

39.4 

1  0.6 

792 

3759.  45 

3596. 86 

27.  1 

1  5.  1 

793 

3808.  16 

3722. 44 

35.  6 

11.6 

794 

3851. 17 

3827.55 

42.7 

9.9 

795 

3879.16 

3942.  93 

49.  6 

8.7 

796 

3896. 63 

4068. 89 

56.  5 

7.9 

797 

3903.40 

4184.  15 

62.  5 

7.3 

798 

3916. 60 

4306.  3  1 

68.8 

6.9 

799 

3904.40 

4428.  30 

74.  5 

6.7 

800 

3879.66 

4554.  38 

80.  2 

6.5 

801 

3851.  19 

4684. 1 3 

85.9 

6.  5 

802 

3810.  16 

4796. 95 

90.7 

6.4 

803 

3752.88 

4908. 69 

95.5 

6.5 

804 

3706. 66 

4998.96 

99.4 

6.  5 

805 

3781. 90 

4945.  89 

97.0 

6.5 

79 


RANGE-  RANG  E  ACCURACIES  (CONTINUED) 


Fix 

X 

Y 

Angle  of 

Radius  of 

No. 

Coordinate 

Coordinate 

Intersection 

90S  Circle 

120 

3261. 07 

411  1.  17 

42.0 

10.0 

12  1 

3262. 5 1 

4368. 76 

63.7 

7.3 

122 

3175.94 

4613. 63 

83. 8 

6.5 

123 

3  189.26 

4667.90 

87.  5 

6.4 

124 

3286.57 

4471.03 

71.9 

6.8 

125 

3317. 76 

4240. 01 

53.9 

8.  1 

126 

3282. 62 

4020. 86 

35.  3 

11.7 

127 

3251.48 

3951.  93 

28.6 

14.3 

128 

3309. 62 

3962. 87 

31.8 

13. C 

129 

3362. 57 

4200. 72 

51.7 

8.4 

130 

3331.06 

4454. 74 

70.7 

6.  9 

13  1 

3252.99 

4654.02 

85.4 

6.5 

132 

3286.69 

4673. 74 

86.3 

6.5 

133 

3366.26 

4474.82 

72.2 

6.8 

134 

3387.38 

4275.  0  1 

57.7 

7.8 

135 

3376.70 

4076.  27 

42.6 

9.9 

136 

3337.33 

3910.52 

29.  0 

14.1 

137 

3345. 86 

3873. 34 

26.7 

15.  3 

138 

3422. 73 

4077. 78 

44.0 

9.6 

139 

3431.43 

4314. 80 

61.3 

7.4 

140 

3396.21 

4530..  82 

76.  0 

6.6 

14  1 

3326. 96 

4700.26 

37.4 

6.  4 

142 

3361.30 

4694. 84 

86.7 

6.4 

143 

3436. 64 

4506.  04 

74 . 4 

6.7 

144 

3474.89 

4309. 23 

61.5 

7.4 

145 

3464. 73 

4120. 80 

48.2 

8.  9 

146 

3425.35 

3933.  68 

34.  0 

12.2 

147 

3459.34 

3928. 23 

34.8 

11.9 

148 

3513.95 

4137.04 

50.  6 

8.5 

149 

3512. 94 

4390. 84 

67.4 

7.0 

150 

3473.67 

4567. 85 

78.4 

6.6 

15  1 

3385. 01 

4766. 90 

90.7 

6.4 

152 

3390.09 

4801.  65 

92. 5 

6.4 

153 

3472. 78 

4657. 67 

83.8 

6.5 

154 

3527.62 

4492. 63 

73.9 

6.7 

155 

3555. 04 

4377.  14 

67.0 

7.  1 

156 

3561. 04 

4288.00 

.61.4 

7.4 

157 

3543. 95 

4204.99 

55.7 

7.9 

158 

3543.69 

4262.  39 

59.5 

7.6 

159 

3543. 85 

4068. 09 

46.8 

9.  1 

160 

3493. 62 

3880.  1  8 

32.9 

12.5 

719 

3519. 28 

3843.  16 

31.6 

13.0 

720 

3586.31 

4061.  28 

47.6 

9.0 

72  1 

3589. 75 

4237. 59 

58.7 

7.7 

722 

3568.78 

4441.  15 

71.1 

6.8 

723 

3525.50 

4623. 58 

81.7 

6.5 

724 

3442. 67 

4796. 98 

91.7 

6.4 

725 

3478.  48 

4818. 98 

92.6 

6.4 

726 

3558.35 

4658.  77 

83.  7 

6.5 

727 

3612.48 

4473.89 

73.4 

6.7 

728 

3635. 37 

4279. 04 

62.  1 

7.4 

729 

3618. 98 

4069.  08 

49.0 

8.8 

730 

3574. 84 

3866.  9  1 

35.  1 

11.8 

73  1 

3587. 72 

3809.47 

32.  1 

12.8 

732 

3613. 95 

3910.  30 

39.  1 

10.7 

733 

3648.48 

4007. 29 

46.1 

9.3 

734 

3676.20 

4121.4  1 

53.  6 

8.2 

735 

3677.  60 

4225.  9  1 

59.8 

7.6 

736 

3676.79 

4332. 54 

66.0 

7.  1 

737 

3658. 87 

4436.  08 

71.7 

6.6 

738 

3638.  19 

4545. 66 

77.6 

6.6 

739 

3604. 38 

4649. 50 

83.2 

6.5 

740 

3572. 48 

4754.  1  8 

88.7 

6.4 

74  1 

3514.26 

4836.98 

93.2 

6.4 

70 


RANGE-RANGE  ACCURACIES  (CONTINUED) 


Fix 

X 

Y 

Angle  of 

Radius  of 

No. 

Coordinate 

Coordinate 

Int  ersection 

90S  Circle 

5  1 

3 1 14. 26 

4776. 58 

96.5 

6.  5 

5? 

2921. 85 

4647.  63 

94. 8 

6.  5 

S3 

2753. 76 

4502.  52 

91 . 5 

6.4 

54 

2793.3  1 

4487. 76 

86.6 

6.  4 

55 

2960. 81 

4619.38 

90.6 

6.4 

56 

3160.  14 

4750.25 

93.6 

6.  4 

57 

3325.  11 

4847.99 

96.  0 

6.5 

58 

3492. 32 

4942.05 

98.4 

6.  5 

59 

3693. 83 

5032. 60 

100.  8 

6.6 

6  1 

3568.  42 

4910.65 

96.3 

6.5 

6? 

3370.43 

48  12. 54 

93.4 

6.4 

63 

3174. 72 

4695.  4  1 

89.  7 

6.4 

64 

2998.35 

4566.  29 

84.5 

6.  5 

65 

2824.36 

443  1.  37 

76.9 

6.6 

6  6 

2819.19 

4368.35 

67.  7 

7.0 

67 

3002. 82 

4519. 53 

80.0 

6.  5 

68 

3184. 74 

4657. 79 

86.  9 

6.  4 

69 

3378. 02 

4763.00 

90.5 

6.  4 

70 

3578. 71 

4871.32 

94.4 

6.5 

71 

3788. 85 

4962. 48 

97.7 

6.  5 

72 

4010. 25 

5031. 39 

100.  4 

6.6 

7  3 

3732. 84 

4895.  17 

95.0 

6.5 

74 

3509.99 

4799.08 

91.3 

6.4 

75 

3307.  18 

4688. 48 

86.9 

6.  4 

76 

3094. 50 

4546.  29 

80.  1 

6.5 

77 

2901.  06 

4391.  68 

68.5 

7.0 

78 

2754.33 

4263. 59 

49.5 

8.7 

79 

3754.  15 

5020.  93 

100.  1 

6.5 

80 

3514.53 

4919.46 

97.2 

6.5 

81 

3287.25 

4788.27 

93.3 

6.4 

82 

3062.52 

4651. 12 

89.4 

6.4 

84 

2846. 46 

4482. 96 

82.4 

6.5 

84 

2845-  13 

4562.17 

91.3 

6.4 

85 

2995. 88 

4403.  69 

68.0 

7.  0 

86 

3165. 76 

4254.99 

52.  7 

8.3 

8  7 

3331.  61 

4  102. 46 

43.3 

9.8 

88 

3500. 54 

3946.  44 

37.  5 

11.1 

89 

3672. 90 

3785.70 

34.0 

12.  2 

90 

3846.42 

3623.  19 

32.0 

12.9 

9  1 

4022. 93 

3464.  18 

31.6 

13.0 

92 

4406. 76 

3998. 00 

68.  3 

7.0 

9  3 

4263.63 

4  159. 42 

70.3 

6.9 

94 

4112. 66 

4331. 20 

73-9 

6.7 

95 

3959.  57 

4527.  77 

79.9 

6.5 

96 

3810. 03 

4631. 47 

85.5 

6.  5 

97 

3672.32 

4864.  18 

93.7 

6.  5 

98 

3535.90 

5040. 98 

102.  3 

6.6 

100 

3065.  44 

4352. 97 

61.3 

7.4 

101 

3068. 28 

4166.32 

41.6 

10.  1 

102 

3077. 08 

4046. 24 

29.5 

13.9 

103 

31 17. 61 

4262. 88 

52.  8 

8.3 

104 

3062.  15 

4499.  1  1 

76.5 

6.6 

105 

3011. 82 

4658.04 

91.7 

6.4 

106 

3122.9  1 

4453.  24 

71.4 

6.8 

107 

3  152.  10 

4199.54 

47.  1 

9.  1 

108 

3118.  73 

4025.  84 

29.3 

14.0 

110 

3188.79 

4154. 61 

43.  8 

9.7 

1  1  1 

3 177.20 

4396.30 

65.9 

7.  1 

112 

3111.79 

4583. 07 

82.8 

6.5 

1  13 

3182. 67 

4534. 87 

77.6 

6.6 

114 

3099.29 

4691.03 

91.3 

6.4 

1  15 

3209.60 

4484.20 

73.3 

6.7 

116 

3235.42 

4271.  54 

55.  2 

8.0 

117 

32C4. 85 

4043. 52 

34.3 

12.0 

APPENDIX  B 

ACCOBACX  CLASSIFICATION :  90  PERCENT  CONFIDENCE  CIRCLES 


CLASSIFIED  RANGE-RANGE  POSITIONS 


Control 

Stations: 

BEACH  LAB  1982 

and  MUSSEL 

1932 

Standard 

Error  Used 

in  Computations:  3  meters 

Fix 

X 

Y 

Anale  of 

Radius 

No.  Coordinate 

Coordinate  Intersection 

90S  £i. 

1 

2668.  05 

4942.  07 

127.  0 

8.2 

2 

2852.  01 

5076.  69 

121.6 

7.7 

•a 

3041. 58 

5194. 85 

113.6 

7.4 

4 

3040. 27 

5148. 74 

117.  0 

7.3 

5 

2838. 26 

5021.35 

120.2 

7.6 

6 

2640.  06 

4877.  13 

126.  9 

8.2 

7 

2553.46 

4752.36 

130.  4 

8.7 

8 

2724.  10 

4885.47 

121.3 

7.7 

9 

2889.06 

5013.  63 

1  17.6 

7.3 

10 

3075.61 

5124.  37 

115.0 

7.2 

1  1 

3 172. 63 

5136.  23 

112.7 

7.0 

12 

2958.  13 

5006.65 

114.  4 

7.  1 

13 

2771. 33 

4876. 02 

117.9 

7.4 

14 

2581. 02 

4729. 89 

126.  9 

8.2 

15 

2584.  4  1 

4665. 31 

124.0 

7.9 

16 

2740. 05 

4805. 63 

1 16.7 

7.3 

17 

29 13. 27 

4928.51 

1 12.8 

7.0 

18 

3097.44 

5047.  77 

111.2 

6.9 

19 

3193.94 

5103.  46 

110.9 

6.9 

20 

2627.47 

4659.  05 

118.8 

7.4 

2  1 

2762.  15 

4769.60 

113.2 

7.  1 

22 

2904. 91 

4875.24 

110.5 

6.9 

23 

3056. 38 

4975. 73 

109.3 

6.9 

24 

3207. 36 

5064.00 

108.9 

6.  8 

25 

3373. 36 

5101.96 

107.  1 

6.8 

26 

3190. 64 

5007. 38 

106.  8 

6.8 

27 

30  15.  45 

4900. 75 

107.  1 

6.8 

28 

2839.  38 

4778.  3  1 

108.  7 

6.8 

29 

2679.99 

4651.  31 

112.8 

7.0 

30 

2727.88 

4632.  77 

106.  9 

6.8 

3  1 

2699. 83 

4776. 97 

105.  2 

6.7 

32 

3092.77 

490  1.  75 

104. 4 

6.6 

33 

3295.42 

5019. 92 

105.0 

6.7 

34 

3502. 07 

5123. 36 

106.  0 

6.7 

35 

3697.02 

5151.09 

105.  3 

6.7 

36 

3474. 77 

5065.  57 

104.  0 

6.  7 

37 

3257.48 

4953.  37 

102.  7 

6.6 

38 

3043.03 

4821. 88 

101.7 

6.6 

39 

2845. 07 

4680.80 

101.7 

6.6 

40 

2746. 96 

4608. 32 

103.  1 

6.6 

4  1 

2748.  09 

4550.40 

97.5 

6.5 

42 

2931.78 

4701.  21 

98.4 

6.5 

43 

3  134.  33 

4834.56 

99.4 

6.5 

44 

3317. 94 

4944. 55 

101.0 

6.6 

45 

35  15.  14 

5043. 56 

102.6 

6.6 

46 

3720.  18 

5117.  15 

103.9 

6.6 

48 

3724. 93 

5083.  62 

102.6 

6.6 

49 

3514.  33 

4989.  13 

100.3 

6.6 

50 

3312. 88 

489  1.  55 

98.5 

6.5 

76 


ELSE  IF 

ELSE  IF 

ELSE  IF 

ELSE  IF 

ELSE  IF 

ELSE 

END  IF 
R  ADIU  S= 
SGX90=2 
SGY90=2 
CIfiAR=3 
EL AR=  3 . 
f ETCJRN 
ENE 


J  *  '  K 


B= 

lic. 

<i£ 


0-710  1*0+1.36546 
.IE. 0.6)  -AND.  (C.GT.C 
C.  544 4*0  +  1 .46488 
.IE.  0. 5)  -  AND.  (C. GT. C 
0.3790*0+1.54/58 
.IE. 0.4)  .AND.  (C. GT. C 
0.2535*0+  1.59778 
.IE. 0.3)  .AND.  (C.GT.C 
0.1652*0+1.62427 
.IE. 0.2) .AND. (C.GT.C 
.094*0+1.63851 


B=  0.0  306*0  +  1 .  64485 


3*SI0X 
. 1 4  6  *  SIG  X 
.146  *SIGY 

.  14 1 5S26*RADI0S**2 
1415926*SGX90*SGY90 


.5) ) THEN 
.4))  THEN 
.3))  THEN 
.2))  THEN 
.1))  THEN 


75 


nn  n  rift  nn  n  nnnnnnnnnnnnn  nnnnnnn 


APPENDIX  A 

SUERCUTINE  FOB  90  EERCENT  CONFIDENCE  CIRCLE  PARAMETERS 


SUERCUTINE  PROE  (SI G 1 . S IG 2 , COR , T 3ET A,  SGXSO,  £G Y  9  0 , 

*  RADIUS, ELAR.C1EAR) 

I MP1ICIT  R  E  AL  *  4  (A-H,0-Z) 

COMPOTES  RADIOS  OF  90%  CONFIDENCE  CIRCLE  (SORT.  METHOD  2) 
THIS  SUBROUTINE  WORKS  FOR  CORRELATED  AND  UNCO F.R  ELAT  ED 
LINES  OF  POSITION. 

INPUT  PARAMETERS: 

SIG 1  AND  SIG2-  STANDARD  ERRORS  OF  TWO  LOP’S 
TBETA  -  ANGLE  OF  INTERSECTION  IN  DEGREES 

(0-180  DEG] 

COB  -  CORRELATION  COEFFICIENT 

(USUALLY  ZERO  EXCEPT  FOR  HYPEFECLIC 
OR  SEXTANT  POSITIONING) 

OUTEUT  PARAMETERS: 

SGX90  AND  SGY90-  SEMI-MAJOR  AND  MINOR  AXES 

OF  90%  ERROR  ELLIPSE 

RADIUS  -  RADIUS  OF  90%  CONFIDENCE  CIRCLE 

ELAH  -  AREA  OF  90%  CONFIDENCE  ELLIPSE 

CIRAR  -  AREA  OF  90%  CONFIDENCE  CIRCLE 


WORK  WITH  AN  ANGLE  IESS  THAN  90  DEGREES 
IF  (TBETA  .  GT.9  0.)  BETA  =  180. -T BETA 
IF  (TBETA.  LE.90  .)  3ETA  =  TBETA 
CHANGE  DEGREES  TO  FADIANS 
RAE=. 017453 2* BETA 

TEA  NS FCR NATION  SIG1  AND  SIG2  TO  CORRELATED 
.  CTANGULAR  SYSTEM 


SIGA  =  SQRT  (  (1./  I  (SIN  (RAD)  ) 

*SIG2*COS(RAD)  ♦SIG2**2)  -SIG2*’*2)  ' 

SIGB=SIG2 

TRANSFORM  COR  RELATION  COEFFICIENT  TO  CORRELATED 
RECTANGULAR  COORDINATE  SYSTEM 

A=  (  (SIG2  +  COS  (RAD)  )  /SIG  1)  +COR 

F=  l/SQRT  (  U2*CCE*SIG2*COS  (RAD)  /SIG  1+  (SIG2/SI 
*  (CCS  (RAD)  ) **2) 

COFAE=A*F 


**2)  )* 
G2**2) 


*  (SIG  1**2*  2.  *COF*SIG 


G  1)  +COR 

♦  COS  (RAD)  /SIG  1+  (SIG2/SIG  1)  **2* 


TRANSFORM  TO  UNCOREELATED  RECTANGULAR 
AA=SQRT (  (SIGA**2*SIGB**2)/2) 

CC  =  SQRT{1-  (4*SIGA**2*SIGB**2*  (1-C0RAB**2)  )/ 
♦  (SIGA**2+SIGB**2)  **2) 

DD=SQET (1+CC) 

SIGX=  AA  *DD 

SIGY  =  SQRT  (SIGA**2+SIGB**2-SIGX**2) 


1* 


COMPUTE  ECCENTRICITY  OF  ELLIPSE 
C=  SIG  Y/SIG  X 

C  COMPUTE  BURT’S  K  FACTOR  BY  LINEAR  INTERPOLATION 
IF  (  (C.  LE.  1.  )  ,AND.(C.GT.  0.9)  )  THEN 
B=  1. 0361  *C>  1.10937 

ELSE  IF  {(C.LE  .  0  .  9)  .  A  N  D.  ( C.  GT  .  0 . 8)  )  TH  S  N 
B  =  0. 9475*C+ 1  .  1896  1 

ELSE  IF  (  (C.IE. 0.8)  -AND.  (C.GT. 0.7) )  THEN 
B=0.8508*C*1 . 26697 

ELSE  IF  ((C.IE. 0.7)  .AND.  (C.GT. 0.6))  THEN 
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Many  variables  exist  when  considering  accuracy  require¬ 
ments  fcr  a  hydrographic  survey.  In  general,  higher  accu¬ 
racy  means  more  time,  money,  and  effort.  Azimuth-a 2imuth 
geometry  is  the  most  accurate  method  of  positioning  analyzed 
in  this  thesis.  This  method  involves  at  least  two  people 
ashore  and  good  ship-to-shore  communications.  Currently, 

NOS  acquires  these  data  manually,  which  minimizes  the  speed 
that  the  vessel  can  operate  and  adds  to  processing  time.  On 
the  otter  hand,  a  survey  using  a  medium-range  system  needs 
little  shore  support  and  the  data  acquisition  is  automated. 
Accuracy  predictions  help  keep  a  balance  between  accuracy 
and  effort.  If  the  desired  accuracy  is  attainable  using  a 
range-range  system  instead  of  an  azimuth-azimuth  system, 
then  the  choice  is  obvious. 

Hydrographic  positioning  in  the  future  will  be  dcmir.ated 
hy  two  methods.  For  cffshore  surveys,  the  Global 
Positioning  System  (GES)  is  expected  to  give  positional 
accuracy  to  10  meters  or  better.  GPS  is  a  satellite  posi¬ 
tioning  system  currently  being  deployed  by  the  Department  of 
Defense  and  will  provide  near  worldwide  coverage  for  users. 
Since  the  full  constellation  of  18  satellites- will  net  he 
operational  until  1988,  it  is  not  yet  known  if  the  expected 
accuracy  of  10  meters  will  be  met.  Nearshore  surveys  may- 
use  multiple  LOP's  for  establishing  hydrographic  positions. 
The  principle  of  least  squares  is  applied  to  redundant 
observations  yielding  the  most  probable  position.  Fcr  both 
GPS  and  least  squares  positioning,  confidence  ellipses  and 
circles  can  be  determined,  although  the  technigues  involved 
are  much  mere  complicated  than  those  presented  in  this 
thesis. 

The  accuracy  classification  scheme  presented  in  this 
thesis  is  predicated  cn  the  elimination  of  systematic 
errors.  Much  work  is  needed  in  identifying  the  sources  of 
systematic  errors  associated  with  hydrographic  positioning 
equipment. 
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lengthened.  In  an  investigation  such  as  this,  it  is  advi¬ 
sable  tc  be  conservative  and  use  the  maximum  length  of  line 
which  is  of era tionall y  feasible  to  provide  coverage  of  an 
area  as  large  as  possible.  The  radius  of  the  90  percent 
confidence  circle  gives  the  hydrographer  a  rough  figure  fcr 
answering  the  guesticn:  Does  the  submerged  pile  exist? 

Knowing  the  parameters  of  the  error  ellipse  could  be 
useful  fcr  conducting  wire-drag,  wire-sweep,  and  side  scan 
sonar  operations.  Fcr  a  position  obtained  with  low  preci¬ 
sion  positioning  equipment,  the  search  to  relocate  a 
-ubmerced  feature  could  cover  a  large  area.  Knowing  the 
parameters  cf  the  error  ellipse  could  reduce  the  area,  time, 
and  effort  of  the  search.  The  search  pattern  could  be 
planned  to  cover  the  desired  confidence  ellipse. 

With  the  quantification  of  accuracy,  a  decision  must  be 
made  concerning  how  much  confidence  is  needed  to  delete  a 
certain  feature  from  the  chart  after  a  search  has  teen  made. 
The  90  percent  confidence  level  may  be  too  low,  whereas  the 
S5  or  99  percent  level  may  suffice.  A  balance  must  be  main¬ 
tained  between  confidence  of  disproval  and  time  and  effort 
spent  on  the  search. 

Accuracy  predictions  in  the  form  of  reliability  contours 
can  be  displayed  using  computer  graphic  terminals.  These 
displays  will  contribute  to  the  efficient  planning  of 
surveys  to  meet  specifications.  Given  the  survey  area,  the 
available  control,  the  positioning  methods,  and  the  preci¬ 
sion  of  the  positionirg  eguipment,  the  hydrographer  can  plan 
the  accuracy  of  the  survey  before  it  is  conducted.  The 
survey  area  and  the  available  control  may  be  such  that  there 
is  flexibility  to  change  control  stations  to  optimize  accu¬ 
racy  ever  an  area  of  critical  importance.  This  information 
can  be  displayed  graphically  and  plans  for  the  survey  can  be 
made  accordingly.  Likewise,  given  an  accuracy  limit,  such 
as  a  10-meter  radius  cf  the  90  percent  confidence  circle, 
the  area  tc  be  covered  at  that  accuracy  can  be  maximized. 
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B.  DSES  FOE  ACCDE AC  I  FIGURES 


NCS  is  currently  developing  the  Shipboard  Data  System 
III  (SDS  III)/  a  hydrographic  data  acguisition  and 
processing  system  which  will  replace  the  present  HYDROLCG/ 
HYDRCILCT  system.  SEE  III  will  revolutionize  data  acguisi¬ 
tion  and  processing  techniques  with  the  capability  tc 
perform  high-speed  calculations  and  display  color  graphics. 
With  this  increased  computer  potential,  data  manipulat icns-- 
such  as  accuracy  computations--can  be  performed. 

Each  position  in  a  survey  can  be  given  a  quality  figure 
based  on  the  radius  of  the  90  percent  confidence  circle. 

This  figure  is  sufficient  for  non-critical  positions  of 
ordinary  hydrographic  data.  Critical  positions  are  these 
which  are  determined  for  significant  features  (i.e.,  wrecks, 
least  depths,  rocks,  and  other  potential  hazards).  For 
these  positions,  the  parameters  of  the  90  percent  error 
ellipse  can  be  computed,  as  well  as  the  radius  of  the  9C 
percent  confidence  circle. 

Many  schemes  can  he  envisioned  for  the  use  of  an  accu¬ 
racy  figure.  For  example,  suppose  the  position  of  a 
submerged  pile  was  determined  by  range-azimuth  geometry  in  a 
prior  survey.  The  radius  of  the  90  percent  confidence 
circle  is  then  5.1  meters  (Ex.  2,  Ch.  II).  The  charting 
agency  now  wishes  to  relocate  the  pile  to  determine  if  it 
still  exists  and  is  still  a  hazard  to  navigation.  In  lew 
water  visibility,  a  common  technique  used  to  resolve  such  an 
item  would  be  to  send  divers  down  over  the  reported  position 
and  conduct  a  circle  search.  Cne  diver  remains  at  the 
reported  position,  holding  a  line,  while  the  other  diver 
swims  a  circumference  holding  the  other  end  of  the  line. 
Theoretically,  if  the  line  is  about  5  meters  long  and  a  hang 
does  ret  occur,  it  is  SO  percent  certain  that  the  pile  has 
been  removed.  For  a  higher  confidence,  the  line  is 
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TABLE  XX 


Accuracy  Pigures  for  o 


Angle  of 
In  ter. 

(deg) 

90 
35 

75 
70 
65 
60 
55 
SO 
95 

30 
25 
20 
15 
10 
5 


K  a 

Ko 

Radius  of 

X 

V 

901  Circ. 

(a) 

(a) 

(a) 

6. 9 

6.9 

6.9 

6.7 

6.2 

.6.5 

7.  1 

6.0 

6.5 

7.5 

5.7 

6.7 

7.9 

5.6 

6.9 

8.  5 

5.  9 

7.2 

9.  1 

§•2 

7.5 

9.9 

5.  1 

8.0 

10.8 

5.0 

8.7 

11.9 

9.9 

9.9 

13.  3 

9.8 

10.  5 

15.  1 

9.8 

11.8 

17.6 

9.7 

13.  7 

2  1.0 

9.7 

16.3 

26.2 

9.6 

20. 2 

39.9 

9.6 

26.8 

52.  2 

9.6 

90.  1 

109.  9 

9.6 

80.  1 

3  a,  o 

12 

Area  of 
Ellipse 

(sg  a) 

130 

131 

i1!! 
139 
1  99 

1  50 
159 
170 
139 

2  60 
309 
391 
503 
750 

1,999 


=  0 


Area  cf 
Circle 

(sg  a) 

130 

131 
135 
191 
199 
162 
179 
203 
235 
280 
395 
990 
588 
832 

1,  289 
2,262 
5,052 
20,  135 


2.  196  for  9GX  probability 


TABLE  X 


Accuracy  Figures  for 

=  10 

*  0 

Angle  of  Ko 

Ko 

Radius  of 

Area  of 

Area  of 

Inter. 

X 

y 

9 OX  Circ. 

Ellipse 

Circle 

(deg) 

(a) 

(a) 

(B) 

(sg  m) 

(sg  a) 

90 

21.5 

21.5 

21.5 

1, 

947 

1,44  7 

85 

22.5 

20.6 

21.6 

1, 

452 

1,459 

80 

23.6 

19.8 

21.8 

1, 

469 

1,497 

75 

29.9 

19.  1 

22.3 

1, 

498 

1^562 

70 

29.5 

18.5 

23.0 

1, 

540 

1,657 

65 

23.2 

18.0 

23.9 

1 

596 

1  798 

60 

30.3 

17.5 

25.2 

l' 

671 

1,989 

55 

32.9 

17.  1 

26.8 

1, 

76  6 

2,252 

50 

35.9 

16.7 

28. 9 

1, 

389 

2,6  15 

95 

39.6 

16.9 

31.5 

2, 

046 

3,117 

90 

99 . 9 

16.  1 

39.9 

2, 

251 

3,335 

35 

50.5 

15.9 

39.5 

2, 

522 

4,890 

30 

58.6 

15.7 

95.  6 

2, 

894 

6,528 

25 

70.1 

15.5 

59.3 

3, 

423 

9,250 

20 

87.9 

15.9 

67.9 

9, 

230 

14,270 

15 

1 16. 3 

15.3 

89.9 

5, 

590 

25,  126 

10 

179.  1 

15.2 

133.7 

8, 

332 

56,  130 

5 

397. 9 

15.2 

266.  9 

16, 

600 

223,727 

R  =  2.196  for  90*  probability 


figures  as  a  function  of  3  for  uncorrelated  TCP's  have  teen 
compiled  using  standard  errors  of  1.3  meters  for  azimuth- 
azimuth  (Table  VIII)  ,  3.0  meters  for  range-range  short-range 
(Table  IX) ,  and  10  meters  for  range-range  medium-range 
(Table  X)  positioning  systems. 


TABLE  Till 

Accuracy  Figures  for  a  -  o  =  1.3  m,  o  =0 

1  2  12 


Angle  cJ 

:  Ko 

Ko 

Radius  of 

Area  of 

Area  of 

Int  er . 

X 

y 

90%  Circ. 

Ellipse 

Circle 

(deg) 

(m) 

(m) 

(m) 

(sg  m) 

(sq  m) 

90 

2.8  I 

2.8 

2.8  | 

24  I 

24 

85 

2.9  i 

2.7 

2.8  | 

1  25 

25 

80 

3.1 

2.6 

2.8  | 

25 

25 

75 

3.2  ] 

2.5 

2.9  i 

!  25 

26 

70 

3.4 

2.4 

3.0  1 

26 

28 

65 

3.7 

2.4 

3.  1  ' 

27 

30 

60 

3.9 

2.  3 

3.3  | 

28 

34 

55 

4.7 

2.3 

3.5  1 

1  30 

38 

50 

5.  1  I 

2.2 

3.8 

i  32 

44 

45 

5.  1 

2.  1 

4.  1 

1  35 

53 

40 

5.8  j 

2.  1 

4.5 

I  38 

65 

35 

6.6 

2.  1 

5.  1 

43 

83 

30 

|  7.6 

2.0  ; 

5.9 

49 

110 

25 

9.1  1 

2.0 

7.1 

58 

156 

20 

11.4 

2.0 

8.8 

71 

241 

1  5 

1  15.1  | 

2.0 

11.6 

94 

425 

10 

S  22.6  I 

2.0 

17.4 

14  1 

949 

5 

45.2 

2.0 

34.7 

28  1 

3,  78  1 

K  =  2.146  for  90!*  probability 


However,  as  mentioned  for  1  :  10,000-scale  surveys  in  a  range- 
range  mode  (a=  3  meters),  this  rule  does  not  always  hold. 

Cn  the  other  hand,  it  is  possible  to  have  g  's  of  less  than 
30°  and  still  meet  specifications.  For  example,  azimuth- 
azimuth  positioning  can  theoretically  be  used  for  B's  of  18° 
to  162°  for  a  1 :  10, 0 CC-scale  survey.  However,  the  eccen¬ 
tricity  cf  the  error  ellipse  is  so  small  that  the  distortion 
introduced  by  using  confidence  circles  can  become 
misleading.  In  view  cf  this,  eccentricities  of  less  that 
0.2  should  not  be  used. 

Using  the  90  percent  radius  criterion,  a  table  has  teen 
assembled  illustrating  the  8  limit  for  various  positioning 
geometries  at  different  survey  scales,  using  assumed  stan¬ 
dard  errors  (Table  VII) .  The  information  in  Table  VII 
illustrates  that  the  20°  to  150°  8  limit  need  not  be  fixed. 
The  B  limits  should  vary  based  on  the  scale  of  the  survey 
and  the  precision  of  the  positioning  equipment.  Accuracy 


Survey 

Scale 


1:2,300 
1:5,000 
1  :  1C, COO 
1 :  2  C , 0  00 
1  : 40, 000 


TABLE  fll 

8  Limits  for  Surveys 


90% 

Radius 


(m) 

2.5 

5.0 

10.0 

20.0 

40.0 


8 


R-R 

1°  =  3) 
Limit 


(deg) 


4  2-138 
27-  153 
23-157* 


R-R 

Jo  t  10) 
8  Limit 

(deg) 


3  5-145 


Az- Az 
(  a  =  1.3) 
8  Limit 

(deg) 


35-145 

23-157* 

23-157* 

23-157* 


*  Eccentricity  limit  of  0.2 

Note:  90%  radii  cf  all  range-azimuth  positions  are 

o=3  and  a 
1  2 


assumed  to  be  5.  1  meters  for 


=  1.3. 
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AZIMUTH-AZIMUTH  ACCURACIES  (CONTINUED) 


X 

Y 

Angle  of 

Radius  cf 

Coordinate 

Coordinate 

Intersection 

90 %  Circle 

4611.20 

4421.29 

91.0 

2.8 

4610. 07 

43  15.  24 

94.  1 

2.  8 

4609. 84 

420  4.  6  6 

97.  4 

2.8 

4603. 99 

4098.63 

100.7 

2.  8 

4600. 52 

3992.  1  b 

104.0 

2.  ? 

4601.27 

3883.  40 

107.3 

2.9 

4601  .  13 

3780.02 

110.  4 

3.0 

4602. 22 

3675.  4  9 

113.  5 

3.  1 

460 1. 48 

3574.  94 

116.  5 

3.  1 

4603.  2  1 

3458. 09 

120.  0 

3.  3 

4524. 84 

3728. 73 

114.8 

3.  1 

4657. 06 

3506.  75 

116.  1 

3.  1 

4648. 60 

3602. 24 

113.  7 

3.  1 

4630.  15 

3696.  1  3 

111.8 

3.0 

4629. 30 

3793.39 

103.  9 

3.0 

4623. 03 

3889.32 

106.  3 

2.9 

4622. 83 

3978.  70 

103.  7 

2.  9 

4617.60 

407  1. 98 

101  .  1 

2.8 

4623.30 

4163.  73 

98.2 

2.8 

4618.  23 

4256.25 

95.7 

2.  8 

82 


CLASSIFIED  RANGE- AZIMUTH  POSITIONS 
Contrcl  Stations:  MISSEL  1932  occupied,  initial  USE  MCN  1978 
Standard  Errors:  Ranee — 3  meters;  T-2 — 1.3  meters 


Fix 

X 

Y 

Angle  of 

Radius  cf 

No. 

Coordinate 

Ccordina  te 

Intersection 

90%  Circle 

414 

2898.58 

4611. 99 

90.0 

5.  1 

4  15 

2979.40 

4483.  3  1 

90.  0 

5.  1 

416 

3027.69 

4342.  47 

90.0 

5.  1 

4  17 

2795.26 

4216. 04 

90.  0 

5.  1 

418 

3002.90 

406  1.  15 

90.0 

5.  1 

419 

2939. 58 

3938. 67 

90.0 

5.  1 

420 

2904. 53 

3888.  38 

90.0 

5.  1 

421 

3195.47 

3380. 99 

90.0 

5.  1 

422 

3128.22 

3487.  82 

90.  0 

5.  1 

423 

3048.  17 

3599.79 

90.0 

5.  1 

424 

2994. 68 

3723. 65 

90.0 

5.  1 

425 

2911. 6  1 

3833.  20 

90.0 

5.  1 

426 

2833.62 

3921.  16 

90.0 

5.  1 

427 

2749.  17 

4000.7  1 

9  0.0 

5.  1 

428 

2752.46 

4029.48 

90.0 

5.  1 

429 

2823. 93 

3969. 92 

90.0 

5.  1 

430 

2894.60 

3920. 66 

90.  0 

5.  1 

431 

2907. 95 

3839. 53 

90.  0 

5.  1 

432 

2961-44 

390  1.  64 

90.0 

5.  1 

433 

3005. 23 

3980.  1  1 

90.0 

5.  1 

434 

3037.78 

4058.0  1 

90.  0 

5.  1 

435 

3065.09 

4134.  5  1 

90.  0 

5.  1 

436 

3074. 88 

4036.  74 

90.0 

5.  1 

437 

3C48.02 

3953.93 

90.0 

5.  1 

438 

2996.76 

3871. 95 

90.  0 

5-  1 

439 

2943.  10 

3809. 76 

90.0 

5.  1 

440 

3128.59 

4049. 87 

90.0 

5.  1 

44  1 

3096. 20 

3978.  78 

90.  0 

5.  1 

442 

3058.84 

3897.24 

90.0 

5.  1 

443 

3009.  6  1 

3825. 86 

90.0 

5.  1 

444 

2914. 71 

3837.97 

90.0 

5.  1 

4  4  S 

30C8.50 

3755.  28 

90.0 

5.  1 

446 

3064.52 

3839. 64 

90-.  0 

5.  1 

448 

3153.77 

4023.  36 

90.0 

5.  1 

449 

3182.38 

4128. 56 

90.  0 

5.  1 

450 

3210.38 

4048. 54 

90.0 

5.  1 

451 

3179. 45 

3970. 85 

90.0 

5.  1 

452 

3146. 52 

3889. 81 

90.  ) 

5.  1 

453 

3099.79 

3820. 34 

90.0 

5.  1 

454 

3047. 75 

3751.  1  8 

90.0 

5.  1 

456 

3037. 99 

3689. 1 9 

90.0 

5.  1 

457 

3121.  90 

3785.62 

90.0 

5.  1 

458 

3186.09 

3883.  49 

90.0 

5.  1 

459 

3234. 50 

4006.  49 

90.0 

5.  1 

460 

3264. 92 

4128. 69 

90.  0 

5.  1 

46  1 

3079. 28 

3672. 85 

90.0 

5.  1 

462 

3159. 97 

3763.  25 

90.0 

5.  1 

463 

3221.69 

3867.47 

90.0 

5.  1 

464 

3273. 65 

3994.  32 

90.  0 

5.  1 

46  5 

3324.  98 

4005.30 

90.0 

5.  1 

466 

3295.53 

3924.  46 

90.0 

5.  1 

467 

3257. 07 

3846.  60 

90.0 

5.  1 

466 

3218. 14 

3770.87 

90.  0 

5.  1 

469 

3172. 68 

3704. 85 

90.0 

5.  1 

470 

3122.00 

3646. 77 

90.0 

5.  1 

47  1 

3097.64 

3628. 62 

90.  0 

5.  1 

472 

3089.39 

3573. 64 

90.  0 

5.  1 

473 

3  118.  69 

3586.  1  4 

90.0 

5.  1 

474 

3210.  37 

3698.  2  1 

90.0 

5.  1 

475 

3262. 22 

3765.90 

90.0 

5.  1 
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RANGE- AZIPUTH  ACCURACIES  (CONTINUED) 


Fix 

X 

Y 

Angle  of 

Radius  of 

NO. 

Coordinate 

Coordinate 

Int  er sect  ion 

90 %  Circle 

476 

3317.  64 

3873.  69 

90.0 

5.  1 

478 

3366. 64 

3900.  37 

90.0 

5.  1 

479 

3335. 51 

3831. 40 

90.0 

5.1 

480 

3303. 12 

3760.21 

90.0 

5.  1 

481 

3258.28 

3696.16 

90.0 

5.  1 

482 

3212. 24 

3631. 36 

90.0 

5.  1 

483 

3165. 13 

3578.90 

90.0 

5.1 

484 

3129. 45 

3550. 77 

90.0 

5.  1 

435 

3142.92 

3499.37 

90.0 

5.  1 

486 

3214. 10 

3575.  56 

90.0 

5.  1 

437 

3280. 84 

3665. 79 

90.0 

5.  1 

488 

3344.  31 

3753. 54 

90.0 

5.  1 

489 

3395. 65 

3847.  40 

90.0 

5.  1 

490 

3420.41 

39 19.  1  1 

90.0 

5.  1 

49  1 

3466. 94 

3949. 20 

90.0 

5.  1 

49  2 

3439.67 

3856.67 

90.0 

5.  1 

493 

3395.38 

3761.44 

90.0 

5.  1 

494 

3335.71 

3669. 35 

90.0 

5.  1 

495 

3275.38 

3583. 62 

90.0 

5.  1 

496 

3197. 78 

3504.  04 

90.0 

5.  1 

497 

3161.00 

3465. 96 

90.0 

5.  1 

APPENDIX  C 

PROGRAM  FOE  90  PIBCENT  CONFIDENCE  ELLIPSE  PASAMETEBS 


PROGRAM  NAME:  ELLIP 


DESCRIPTION: 


COMPUTES  ORIENTATION ,  90%  SIGHA-X,  90%  SIGMA-Y, 
FOR  ERROR  ELLIPSE  ABOUT  A  HYDROGRAPHIC 
POSITION  ESTABLISHED  BY  RANGE-RANGE,  AZIMUTH- 
AZIMUTH  OS  RANGE-AZIMUTH  POSITION 


AUTHOR:  NICHOLAS  E.  PERUGINI 
LT.  NO AA 

NAVAL  POSTGRADUATE  SCHOOL 
DATE:  SEPTEMBER  ,  1984 


IMPLICIT  REAL  *  4  (A-H,C-Z) 

INITIALIZE  VALUES:  FOR  RANGE-RANGE  AND 
-XL  AND  YL  ARE  COORDINATES 
-XR  AND  YR  ARE  COORDINATES 
-SIGI  AND  SI GR  ARE  RESPECTIVE 
ASSOCIATED  WITH  EACH  LOP 


AZIMUTH-AZIMUTH: 

OF  LEFT  STATION 
OF  RIGHT  STATION 
STANDARD  E-.  ROSS 


C 

C 

c 


FOR  RANGE-AZIMUTH: 

-XL  ANE  YL  ARE  COORDINATES  OF  OCCUPIED  STATION 
-SIGI  IS  SIGMA  OF  THEODOLITE  LOP 
—  S IGR  IS  SIGMA  OF  RANGE  LOP 
************************************************************ 
XL=4914. 75 
YI=2  C09. 86 
SI GL=  3- 0 

XR=2474. 75 
YR=4247. 42 
SIGR=  3-0 

ENTER  CORRELATION  COEFFICIENT:  USUALLY  ZERO  FOR  R-E, 

R-AZ,  AND  AZ-AZ 

RC  =  0.0 
PI=3. 141593 


ENTER  INDICATOR  TO 


TELL 
IND 
IND  = 
IND 


WHAT 
=  1 
2 


KIND  OF  DATA 
RANGE-RANGE 
AZIMUTH-AZIMUTH 
RANGE-AZIMUTH 


IS  ENTERING  PROGRAM 


INC  I 

S  A 

TOGGLE 

WH 

ICH  C 

NOTH! 

RG 

IN  PROGR 

AM 

SHCU 

READ 

IN 

DATA  FRO 

M 

£  ATA 

PX 

EY 

TD 

SENTI 

NEL 

IS  I FIX 

= 

999, 

10 

CONTINUE 

EE 


FOR  BETA  GREATER  THAN 
CHANGED  FROM  HERE  ON 


=  3 

*******************44***************************** 
IND=  1 

C 
C 
C 

c 
c 
c 
c 
c 
c 


90  DEG. 


FILE:  IF IX  =  FIX  NUMBER 
=  X  COORDINATE  OF  HYDRO  POSITION 
=  Y  COORDINATE  OF  HYDRO  POSITION 
=  ANGLE  OF  INTERSECTION  IN  DEGREES 
TELLS  PROGRAM  TO  STOP  READING 


85 


nftnn  onnn  no  n  nnnrifinnn  no 


READ  (4,  20)  IFIX,ILL,PX, FY.TD.RAD 
20  FORMAT (1 X,I3 , 3X, I 1 , oX , r 7. 2,61,71.2,  5X,F9.4,3X,F5.2) 

DC  =  0 

IF  (IFIX. EC. 999)  GO  TO  900 
IF  (TD.LT.90.)  ED=TD 
IFJTD.LT.  90. }  GO  TO  30 

C  WORK  WITH  BETA  LESS  THAN  90  DEGREES:  TOGGLE  TURNED  CN  TC  ONE 
DE=180.-TD 
INC=  1 

30  CONTINUE 
C 

C  KEEP  TANGENT  FUNCTION  FROM  GOING  UNDEFINED  IN  A  RARE  CASE 
C  OF  TEE  FIX  AND  CONTROL  STATION  HAVING  SAME  COORDINATES 


IF  (EX. EQ. XL)  P  X=  PX  +  0. ! 
IF(PY.EQ.YL)  PY=PY+  0.! 


CHANGE  DEGREES  TO  RADIANS 
3ETA  =  . 0174532  9*DD 

USE  LEFT  STATION  AS  BASIS  FOR  COMPUTATIONS 

ORIENTATION  ANGLES  HILL  BE  FIXED  WITH  RESPECT  TO  LEFT  LOP 

FINE  AZIMUTH  FROM  NORTH  BETWEEN  HYDRO  POSITION  AND  IEFT 
STATION.  AZIMUTH  HILL  EE  DEFINED  EETWEEN  0-130  DEGREES 
MEASURED  CLOCKWISE  FROM  NORTH. 

THIS  IS  THE  RANGE-RANGE  AZIMUTH  DETERMINATION. 

IF  (IND.NE.  1)  GO  TO  40 
IF(PY.GE.YL)  TEEN 

IF(PX.GE.XL)  THEN 

ALPHA  =  PI-ATAN { (PY-YL) / (PX-XL) ) 


EISE 


ELSE 

ALPHA  = 
END  IF 


AT  AN  (  (PY-YL)  /  (XL-PX)  ) 


IF(PX.GE.XL)  THEN 

ALPHA  =  ATAN  (  (YL-PY)  /  (PX-XL)  ) 

ELSE 

ALPHA  =  PI— ATAN ( (YL-PY)/(XL-PX)) 

END  IF 
END  IF 
GO  TO  60 

AZIMUTH  FIXING  FOR  AZIMUTH-AZIMUTH  POSITIONS 
40  CONTINUE 

IF  (EY.GE. YL)  TEEN 

IF(PX.GE.XL)  THEN 

ALPHA  =  ATAN  ( (PX-XL) / (PY- YL) ) 

ELSE 

ALPHA  =  i  I-ATAN  (  (XL-PX)  /  (PY-YL)  ) 

END  IF 

ELSE 

IF(PX.GE.XL)  THEN 

ALPHA  =  PI-ATAN(  (PX-XL)/ (YL-PY)  ) 

ELSE 

ALPHA  =  ATAN { (XL-PX) / (YL-PY) ) 

END  IF 
END  IF 

AZIMUTH  EQUALS  THETA  FOR  RANGE  AZIMUTH  CASE,  ASSUMING 
THEODOLITE  SIGMA  IS  LESS  THAN  RANGE  SIGMA 

IF  (IND.EQ.3)  GO  TO  70 

BEGIN  COMPUTING  THETA,  THAT  IS  THE  ANGLE  OF  ROTATION  FROM 
LEFT  LCP 
60  CONTINUE 

E  1=SIGL**2*SIN  (2*BETA)  +2*RO*SIGL*SIGR*SIN (BETA) 
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nnn  no  nnnnnn 


B2=SIGL**2*COS  |2*BE?A)  +2*R0* SIG L*SIGR*COS  (BETA)  +SIGR**2 
IF  e(ABS^B2)  .IT. 0.000  1)  E2=.0001 

C  COMPUTE  ROTATION  ANGIE  FROM  IEFT  LOP 
TH=0.5*ATAN (B3) 

90  CONTINUE 


DEFINE  SEMI-MAJOR  AXES  ORIENTATION  IN  TERMS  OF  0-180  DEGREES 
ROTATION,  CLOCKWISE  FROM  NORTH 

RANGE-RANGE  CASE 

IF(IND.EQ.I)  TEEN 

IF  (INC.  EQ  .  1)  THETA  =  A  IPHA  +  TH 
IF  'INC. EQ. 0)  THETA=AIPHA-TH 
END  IF 

AZIMUTH-AZIMUTH  CASE 

IF  (IND.  EQ.  2)  TEEN 

IF  (INC . EQ . C)  THETA  =  AIPHA  +  TH 
IF(INC.EQ.I)  THETA  =  AIPHA-TH 
END  IF 


FIX  ROTATION  ANGLE  FROM  0-180  DGREES 
70  CONTINUE 

C  CONDITION  FOR  RANGE- AZIMUT H  DATA 
IE  IND. EQ . 3)  THETA=AL?HA 
IF  THETA. LT.0.)  THETA=  EI+THETA 
IF  THETA. GT. PI)  THETA=  THETA-PI 
C  DEG  IS  THE  SEMI-MAJOR  ELLIPSE  AXIS  ORIENTATION 
DEG=  57.  295779  +  THETA 

COMFUTE  90%  SIGMAX  AND  SIGMAY  OF  ERROR  ELLIPSE 


IN  DEGREES 


100 


CALL  PROS  (SIGI.SIGR.RO,  ID,  SGX90.SGY  90.  RADI  US,  EL  A  E,  CIRAR) 

WRIT  E (7 ,  100}  IFIX,?X.PY,TD.SGX50.SGY§0,DEG 

FORMAT flX,I3,3X,F7.2,3X,F7.2,3X,F5.  1 , 3X,  F 4. 1 , 3X , F4 . 1 , 


900 


♦  3X,F5.  t) 
GO  TO  ID 
CONTINUE 
STOP 
END 


APPENDIX  D 

ACCURACY  CLASSIFICATION:  SO  PERCENT  CONFIDENCE  ELLIPSES 


CL ASSI TIED  RANGE-RANGE  POSITIONS 
Control  Stations:  BEACH  LAB  1982  and  MUSSEL  1932 
Standard  Error  Used  ic  Computations:  3  meters 


Fix 

X 

Y 

Beta 

90% 

90% 

Orien  ta- 

No. 

Coord. 

Coord. 

Siqma  X 

Si<jma  Y 

tion 

1 

2668.05 

4942.07 

127.0 

10.2 

5.  1 

79.0 

2 

2852. 01 

5076.69 

121.6 

9.3 

5.2 

85.3 

3 

3041.58 

5194.85 

118.6 

8.9 

5.  3 

90.2 

4 

3040.27 

5148.74 

117.0 

8.7 

5.3 

90.6 

5 

2838.26 

5021.35 

120.3 

9.  1 

5.2 

85.3 

6 

2640.06 

4877.13 

126.9 

10.2 

5.  1 

78.1 

7 

2553.46 

4752.36 

130.4 

10.9 

5.0 

74.1 

8 

2724.  10 

4885.47 

121.4 

9.3 

5.2 

82.0 

9 

2889.06 

5013.63 

1  17.  6 

8.  8 

5.3 

87.2 

10 

3C75.61 

5124.37 

115.0 

8.5 

5.4 

91.9 

1  1 

3172.63 

5136.23 

1  12.7 

8.2 

5.5 

94.5 

12 

2958.  13 

5006.65 

114.4 

8.4 

5.4 

89.7 

13 

2771.33 

4876.02 

1  18.0 

8.8 

5.3 

84.2 

14 

2581.02 

472S.89 

127.0 

10.2 

5.  1 

75.9 

15 

2584.41 

4665.31 

124.0 

9.7 

5.2 

76.7 

16 

2740. 05 

4805.63 

1  16.7 

8.7 

5.3 

83.8 

17 

2913.27 

4928.51 

1  12.8 

8.2 

5.  5 

89.2 

18 

3C97.44 

5047.77 

111.2 

8.  1 

5.5 

93.5 

19 

3193.94 

5  103.46 

110.9 

8.0 

5.5 

95.5 

20 

2627.47 

4659.05 

1  18.8 

8.9 

5.3 

79.8 

21 

2762.15 

4769.60 

113.2 

8.3 

5.5 

85.4 

22 

2904.91 

4875.24 

110.5 

8.0 

5.5 

89.7 

23 

3056.38 

4975.73 

109.3 

7.9 

5.6 

93.3 

24 

3207.36 

5064.00 

108.9 

7.8 

5.6 

96.3 

25 

3373.36 

510 1.96 

107.  1 

7.7 

5.7 

100.0 

26 

3  19  0.64 

5007.38 

1  06.8 

7.6 

5.7 

96.7 

27 

3015.45 

4900.75 

107.1 

7.7 

5.7 

93.2 

28 

2839.38 

4778.31 

108.7 

7.8 

5.6 

88.9 

29 

2679.99 

4651.31 

1  12.8 

8.2 

5.5 

83.4 

30 

2727.88 

4632.77 

1  06.9 

7.6 

5.7 

86.7 

31 

2899.83 

4776.97 

105.2 

7.5 

5.7 

91.3 

32 

3092. 77 

490  1.75 

104.4 

7.4 

5.8 

95.6 

33 

3295.42 

5019.92 

105.  0 

7.5 

5.7 

99.2 

34 

3502.  07 

5123.36 

106.0 

7.6 

5.7 

102.6 

35 

3697.02 

515  1.09 

1  05.3 

7.5 

5.  7 

106.2 

36 

3  47  4.  77 

5065.57 

104.1 

7.4 

5.8 

102.7 

37 

3257.48 

4953.37 

1  02.  7 

7.3 

5.  3 

99.3 

38 

3043.03 

482  1.88 

101.7 

7.  2 

5.9 

95.5 

39 

2845.07 

4680.80 

101.7 

7.2 

5.9 

91.4 

40 

2746.96 

4608.32 

103.1 

7.3 

5.8 

88.6 
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CLASSIFIEE  AZIMUTH-AZ IMUTH  POSITIONS 


Control  Stations:  USE  MON  1978  and  MUSSEL  1932 
Standard  Error  Used  in  Computations:  1.3  meters 


Fix 

X 

Y 

Beta 

9  0^ 

90% 

Orienta¬ 

No. 

Ccord. 

Coord. 

Si^ma  X 

Si^ma  Y 

tion 

619 

4449.26 

271  1.40 

156.9 

9.8 

2.0 

139.4 

620 

4437.86 

2807.29 

153.0 

8.5 

2.0 

139.8 

621 

4427.41 

2897.17 

1  49.6 

7.5 

2.0 

139.  8 

622 

4419.65 

2989.36 

146.2 

6.3  ■ 

2.  1 

129.  8 

623 

441C.42 

3072.20 

143.4 

6.3 

2.  1 

139.  6 

624 

4290.60 

.  3154.04 

141.3 

5.9 

2.  1 

139.1 

625 

4376.52 

3234.24 

138.9 

5.6 

2.  1 

138.  6 

626 

4364.31 

3319.82 

136.2 

5.3 

2.  1 

138.  C 

627 

4348.71 

34  1  1.24 

1  33.5 

5.0 

2.  1 

137.3 

628 

4334.31 

3502.51 

1  30.  7 

4.7 

2.2 

136.5 

63  1 

4238.97 

338  1.55 

135.1 

5.2 

2.  1 

137.4 

632 

4337.95 

3290.68 

138.7 

5.6 

2.  1 

137.  8 

633 

4227.70 

3199.92 

142.8 

6.2 

2.  1 

133.  1 

634 

4322.53 

3107.83 

1  46.9 

6.9 

2.  1 

138.2 

635 

4223.18 

3012.49 

151.0 

7.9 

2.0 

138.  3 

636 

4324.42 

2916.00 

155.3 

9.2 

2.  0 

138.  1 

637 

4327.25 

2818.39 

1  59.8 

11.3 

2.0 

137.7 

638 

4394.77 

2806.86 

156.0 

9.  5 

2.0 

138.9 

639 

4386.58 

2903.09 

1  52.0 

8.2 

2.0 

139.  1 

640 

4377.29 

2998.63 

1  48.4 

7.2 

2.  1 

139.  1 

641 

4367.00 

3090.26 

145.  1 

6.6 

2.1 

138.9 

642 

4355.55 

3187.06 

141.9 

6.  0 

2.  1 

139.5 

643 

4345.97 

3285.88 

138.5 

5.6 

2.1 

138.0 

644 

4256.90 

3516.02 

133.6 

5.0 

2.  1 

1  -C.  c 

645 

4260.68 

3416.65 

137.  4 

5.4 

2.1 

136.  2 

6  46 

4264.77 

332  1.01 

141.  1 

5.9 

2.1 

136.8 

647 

4283.59 

3208.82 

144.8 

6.5 

2.1 

137.5 

648 

4293.24 

3127.36 

147.7 

7.  1 

2.  1 

127.8 

649 

4200.30 

3024.70 

151.9 

8.1 

2.0 

137.9 

650 

4345.53 

3145.58 

144.  1 

6.  4 

2.  1 

138.5 

651 

4370.38 

3236.09 

139.  1 

5.7 

2.  1 

138.  5 

652 

4398.77 

3327.77 

134.2 

5.  1 

2.  1 

128.  4 

653 

441  1.48 

342  1.46 

130.2 

4.7 

2.2 

138.0 

654 

4438.30 

3506.23 

125.9 

4.3 

2.2 

137.7 

655 

4470. 97 

359  1.21 

121.6 

4.0 

2.3 

137.  4 

656 

4502.73 

3677.50 

117.4 

3.8 

2.3 

137. C 

657 

4514.38 

3767.26 

114.0 

3.6 

2.4 

136 . 2 

658 

4512.00 

3860.13 

111.1 

3.  5 

2.4 

135.2 

659 

4520.  10 

3948.53 

107.9 

3.4 

2.4 

134.3 

660 

4494.79 

3049.75 

139.2 

5.7 

2.  1 

141.  C 

CLASSIFIED  RANGE-AZIMUTH  POSITIONS 


Control 

.  Stations: 

0SE  MON 

1  978, 

MUSSEL 

1932 

Standard  Errors: 

RANGE — 3 

me  ters ; 

T-  2 —  1 

.3 

Fi  x 

X 

7 

3eta 

90% 

90% 

Orienta¬ 

No. 

Ccord. 

Coord. 

Si^ma  X 

Siqma  I 

tion 

414 

2E38.53 

4611.99 

90.0 

6.  4 

2.  8 

49.3 

415 

2979.40 

4483.31 

90.0 

6.4 

2.8 

64 .  S 

416 

3027.69 

4342.47 

90.0 

6.4 

2.8 

80.2 

417 

2795.26 

4216.04 

90.0 

6.4 

2.8 

9  5.6 

418 

3CC2.S0 

406  1.15 

90.0 

6.  4 

2.8 

1C9.4 

4  19 

2939.58 

3938.67 

90.0 

6.4 

2.8 

123.  6 

420 

2904.53 

3888.38 

90.0 

6.  4 

2.8 

129.5 

421 

3  195.47 

3380.99 

90.0 

6.4 

2.8 

140.  2 

422 

3128.22 

3437.82 

90.0 

6.4 

2.3 

139.3 

423 

3048.  17 

359  S. 79 

90.0 

6.4 

2.8 

138.5 

424 

2S94.68 

3723.65 

90.0 

6.  4 

2.8 

135.2 

425 

2911.61 

3833.20 

90.0 

6.4 

2.8 

133.  5 

426 

2833.62 

3921.16 

90.0 

6.  4 

2.8 

132.3 

427 

2749.  17 

400C.71 

90.0 

6.4 

2.8 

132.0 

428 

2752.46 

4  029.48 

90.0 

6.  4 

2.  8 

128.1 

429 

2823.93 

3965.92 

90.0 

6.4 

2.8 

128.  5 

430 

2894.60 

3920.66 

90.0 

6.  4 

2.3 

127.9 

43  1 

2507.95 

3839.53 

90.0 

6.  4 

2.8 

133.  3 

43  2 

256  1.44 

390  1.64 

90.0 

6.  4 

2.  8 

125.4 

433 

3C05.23 

398C. 11 

90.0 

6.4 

2.8 

1  16.7 

434 

3C37.73 

4058.01 

90.0 

6.  4 

2.8 

103.6 

435 

3C65.09 

4134.51 

90.0 

6.4 

2.8 

100.  8 

436 

3C74.88 

4036.74 

90.0 

6.4 

2.8 

ICS.  3 

437 

3C48.02 

3953.93 

90.0 

6.4 

2.8 

117.  1 

438 

2996.76 

3871.95 

90.0 

6.  4 

2.8 

125.7 

439 

2943.  10 

380  5.76 

90.0 

6.4 

2.8 

133.  1 

440 

3  128.59 

4049.87 

90.0 

6.  4 

2.8 

106.8 

44  1 

3CS6.20 

3978.78 

90.0 

6.  4 

2.8 
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